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ABSTRACT 


Data from a line of magnetometers stretching Arr 
a corrected geomagnetic meridian ~ 302°E through western 
Canada are used to study the relationship between the convec- 
tion westward eliectrojet and Pc 5 micropulsations inthe 
morning sector. It was found that the dominant spectral 
bands in Pc 5 range occur within the same latitudinal range 
occupied by che electrojet. The intensity contours and the 
character of the polarization parameters clearly show that the 
pulsational activity follows the westward electrojet. 
Evidence of spatial oscillations of the borders of the elec- 
trojet at dominant Pc 5 frequencies is presented. The 
response of the Pc 5 micropulsations to rapid reconfigurations 
of the electrojet and the impulsive nature of the Pc 5 wave 
trains suggest that these micropulsations represent a tran- 
Sient response to changes in the three dimensional magneto- 
Sphere-ionosphere current system associated with large scale 
convection within the magnetosphere. The displacement 
current in the F-region and induction in the earth are 
suggested to have an important effect in the determination 
OTetherpolanizavion of Pc 5) micropulsations: 

In addition, an analysis and interpretation of 
data for a few rare ‘Pc 4* giant-micropulsations in the 


morning sector recorded during 1971 and 1974 (when three of 
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the stations were distributed along the same latitude 

line) indicate that these pulsations are highly localized 
both in latitudinal and longitudinal extent and occur 
during periods of extremely low activity in the magneto- 
sphere. The 'Pc 4' activity appears to be influenced by 
the initiation of substorm activity; however, the region of 
pulsation disturbance is found to lie outside of the elec- 
trojet regions. It is suggested that the 'Pc 4' giant 


pulsations may be due to field line resonances. 
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] Three dimensional representation of the 9 
configuration of the magnetosphere [after 
Heikkila, 1973]. The diagram has been 
modified by Rostoker [1976] to show the 
plasma mantle particle distribution. 


2 Schematic diagram of the magnetosphere 1] 
projected on the noon-midnight meridian 
plane illustrating the Birkeland current 
flow in the evening sector [after Rostoker 
and Bostrém, 1976]. The diagram has been 
modified to show the locations of the 
auroral oval. 


3 An artist's conception of the auroral oval. 12 


4 Electric field configuration in the magne- 132 
. totail projected on the cross-sectional 

area of the tail [after Rostoker and 

Bostrom, 1976]. 


5 Direction of ionospheric currents (Hall 15 
current I, and Pedersen current 1) and 


Birkeland current flow and the orientation 
of the ionospheric electric field in the 
polar ionosphere [after Rostoker and 
Bostrom, 1976]. The diagram has been 
modified to show the flow of currents in 
the dayside and the net downward field- 
aligned currents in the morning sector and 
the net upward field-aligned currents in 
the pre-midnight quadrant. 


6 An example of Pc 5 micropulsations recorded 17 
at the station of Fort Reliance (70.3°N 
geomagnetic). The disturbance is recorded 
in the local magnetic coordinate system 
CHD 2. HH pointse COWdrds, tmecshoca | 
magnetic north, D east and Z downward. 

Local time is approximately UT less seven 
hours . 
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Polar plots of the auroral oval (between 

dotted lines), the intensity maxima of 
Pc 5's (between solid lines), and the 

plasmapause (dashed lines) in corrected 
sf las coordinates [after Samson, 
97 2a)" 


The diurnal and latitudinal variations of 
the; sense? of polarization for’ Pc 5’ pulsa- 
tions of ~ 5 mHz in the horizontal plane 
[after Samson et al., 1971]. 


Average sense of polarization in the 
horizontal plane of micropulsations in 

the spectral range 3.3 - 100 mHz. The 
Shaded regions marked C have predominantly 
clockwise polarization and those marked CC 
have predominantly counter clockwise 
polarization. The unshaded regions have 

no predominant sense of polarization [after 
Rankin and Kurtz, 1970]. 


An example of 'Pc 4' giant pulsations 
recorded at station of Fort McMurray 
(64.2°N geomagnetic). 


Pot are plots’ of the velocity vectors for 
the three waves which can be excited in 
a highly conducting compressible hydro- 
magnetic medium [after Ferraro and 
Plumpton, 1966]. 


Schematic diagram of equatorial plane of 
the magnetosphere showing the senses of 
polarization induced by a Kelvin-Helmholtz 
instability at the magnetopause [after 
Orr, 1973]. 


Projection of the particle trajectory 
showing the diffusions of protons due 

to waves with a magnetic field component 
in) the-direction of the’ principal. nornal 
to the dipole field lines [after Dungey, 
1965]. 


Schematic representation of the variation 
of amplitude and polarization as a func- 
tion of L for unstable surface waves on 
the magnetopause coupled to a resonant 


field line [after Chen and Hasegawa, 1974a]. 
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Latitude-dependence of magnetic variations 65 
during four time periods when sinusoidal 
oscillations were observed at the lowest 
latitude. All three stations are located 

at the same magnetic longitude [after 

Lanzerotti et al., 1973]. 


The local time distribution of the detached 7\ 
plasma regions [after Chappell, 1974]. 


Locations of the magnetometer stations in 74 
corrected geomagnetic coordinates. The 

data recorded by stations indicated by 

Squares were used in this thesis. Data 

from the stations indicated by solid circles 
were used to determine the region over which 

the micropulsation disturbance extended. 

Only standard magnetograms were available 

from the stations indicated by solid cir- 

cles. 


Magnetograms from the University of Alberta 83 
magnetometer stations for the event 

eceursgingson Day 235, 1971. The time ‘sicale 

is marked in hours of Universal Time. The 
station code names and coordinates are 

listed in Table 2. 


1-20 mHz bandpass filtered micropulsation 84 
data sforethes Day 2365 ,0197.) event. .The 
Original data are shown in Figure 18. 


1-20 mHz bandpass filtered micropulsation 85 
data for the Day 235, 1971 event. 

Magnetogram for the Day 16, 1972 event. 86 
1-20 mHz bandpass filtered micropulsation 87 
dataaror thesDay 16541972 event. 

Magnetograms for the Day 17, 19/72 event. 88 
1-20 mHz bandpass filtered micropulsation 89 


datattor the Day 17, 1972 event. 
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Latitude profiles showing the magnetic 
perturbation pattern along the Alberta 
meridian line. The data are presented 
in geomagnetic dipole coordinates with 
X positive northwards, Y positive 
eastwards and Z positive downwards. 
The arrowheads at the bottom of each 
frame delineate the borders of the 
convection westward electrojet. The 
Universal Time at which the profile is 
taken is given at the top of the frame. 


Plots of contours of equal relative 
intensity of micropulsations in the 
1.3-2.5 mHz spectral band (upper panel). 
The location of the electrojet (indicated 
by solid lines) is superposed on the 
intensity maximum indicated by the dark 
band (lower panel) for the Day 16, 1972 
event. The locations of our stations 

are indicated by the stars on the 
Ordinate axis. 


Same as Figure 26 except for the Day 235, 
1971 event. 


Same as Figure 26 except for the Day vhiv, 
1972 event. 


Power spectra of the three magnetic com- 
pone mts aC sD3Z)j% at ithe Uni vens pty sof 
Alberta magnetometer stations. The top 
panel shows the power spectra at our 
most northern station, with the lowest 
panel showing the spectra at our most 
southern station. The spectra were 
calculated in the interval indicated at 
the top of the Figure. The records were 
filtered with a 1-20 mHz bandpass filter. 
The stations that were within the con- 
fines of the westward electrojet during 
most of the time interval are indicated 
by the term WJ. These power spectra 

are for the Day 17, 1972 (1100-1200 UT) 
event. 
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Same as Figure 29 except for the Day 17, 
1972 (1200-1300 UT) event. 


Same as Figure 29 except for the Day 16, 
1972 (1700-1800 UT) event. 


Same as Figure 29 except for the Day 17, 
1972 (1700-1800 UT) event. 


Latitude profiles of the powers in three 
components in different spectral bands. 
The borders of the westward electrojet 
are indicated by the dashed lines. This 
latitude profile is for the Day 16, 1972 
(1500-1600 UT) event. 


Same as Figure 33 except for the Day 17, 
1972 (1500-1600 UT) event. 


Same as Figure 33 except for the Day 17, 
1972 (1700-1800 UT) event. 


Polarization ellipses in the H-D plane. 
The centres of the ellipses are at the 
latitudes of the stations which are 
indicated by the stars on the ordinate 
axis. The sizes of the ellipses are not 
scaled to the average intensities and all 
the major axes have the same length. 
Ellipses that are shaded depict clockwise 
(CW) polarization while the ellipses that 
are not shaded depict counter-clockwise 
(CC) polarization. The location of the 
convection westward electrojet indicated 
by solid lines is superposed in the 
batckignotinids. The. plot as fomither Day! alrer, 
1972 event. 


Same as Figure 36 except that the polari- 
zation ellipses are in the H-Z plane for 
the Day 16, 1972 event. 


Same as Figure 36 except that the polari- 
zation ellipses are in the D-Z plane for 
thewvay lo, 1972 event. 


Same as Figure 36 except for the Day 235, 
1971 event. 7 
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Same as Figure 37 except for the Day 235, 
1971 event. | 


Same as Figure 38 except for the Day 235, 
1971 event. 


Same as Figure 36 except for the Day 17, 
1972 event. 


Same as Figure 37 except for the Day 17, 
1972 event. 


Same as Figure 38 except for the Day 17, 
1972 event. 


Variation in the sense of polarization in 
H-D plane as a function of invariant 
latitude and Universal Time. The loca- 
tion of the electrojet is superposed in 
the background. The locations of our 
stations are indicated by the stars on 
the ordinate axis. CW refers to clock- 
Wise; DOharizatvon. [hs =plot as sstorethe 
Day 17, 1972 event. 


Same as Figure 45 except for polarization 
in the H-Z plane. 


Same as Figure 45 except for polarization 
in the D-Z plane. 


Latitude profile of powers and relative 
phase changes for the Day 17, 1972 
(1500-1600 UT) event at 4.0 mHz. The 
locations of our stations are indicated 
atcthe stop. tINLAT stands for invariant 
latitude. The borders of the convection 
westward electrojet are indicated by the 
dashed lines. 


Same as Figure 48 except for the Day 16, 
1972 (1700-1800 UT) event at 2.2 mHz. 
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59 


Plot of demarcation lines (dotted lines) 
in the H-D plane and the line of intensity 
maximum (dashed line) as a function of 
invariant latitude and Universal Time. 

The location of the convection westward 


electrojet is superposed in the background. 


Ehitsi (prot ris “for the Day 235; -P971 tevent 
in the 2-3 mHz band. 


Same as Figure 50 except for the Day 16, 
1972 event in the 1.3-2.5 mHz band. 


Same as Figure 50 except for the Day 17, 
1972 event in the 3.3-5.1 mHz band. 


Plots of intensity contours and the con- 
vection westward electrojet. The relative 
intensity at each station was determined 
with respect to the maximum intensity for 
the entire interval. These plots are 

for the Day 235, 1971 event in the 2.0- 
3.0 mHz band. 


Same as Figure 53,except for ithe Day 'h6, 
1972 event in the |..3-2.5 mHz band. 


Same as Figure 53 except for the Day 16, 
1972 event in the 4.1-5.3 mHz band. 


Same as Figure 53 except for the Day 17, 
1972 event in the 1.5-2.8 mHz band. 


Same as Figure 53 except for the Day 17, 
1972 event in the 3.3-5.1 mHz band. 


Location of the convection westward elec- 
trojet and the filtered magnetograms for 
CONT and FTCH in the Z-component. The 
borders of the electrojet were determined 
at intervals of 23 sec. For display 
purposes the magnetograms are placed 
between the electrojet borders. This 
plein s.afor athe: Day 17, 1972,{( 1300-1400 
UT) event. 


Same as Figure 58 except that the pole- 
ward border is not shown and the magneto- 
grams are for MCMU and MENK in the H- 
component, fhe plot is for=the Day 17; 
1972 (1200-1400 UT) event. 
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Same as Figure 58 except that the magne- 
tograms are for CONT and RELI in the H- 
components “The plot ts for the Day.17; 
1972 (1700-1800 UT) event. 


Plot showing the oscillations of the 
equatorward border of the convection 
westward electrojet and the magnetogram 
for SMIT for the Day 17, 1972 (1500-1800 
UT) event. The lower panel shows the 
variation in magnitude of the negative 
AX peak in the latitude profiles of 
magnetic perturbations. The units for 
HMAX are y. 


An equivalent LCR circuit for the proposed 
model three-dimensional current system. 

e is the energy source associated with 

= aes of energy from convecting 
plasma. is the inductance associated 
with the nllonetic energy in the generator 
region. is the capacitance associated 
with the ne energy of the convecting 
plasma in the generator region. Cy is 


the capacitance associated with the F- 
region of the ionosphere. Ry is the 


resistance associated with the E-region 
of the ionosphere. 


The diurnal variation of the ionospheric 
electric field in the region of the 
auroral oval [after Mozer and Lucht, 
1974]. 


Field lines from rings of constant 
latitude as viewed from above and pro- 
jected into the equatorial plane. Lines 
are shown for every 2 hours of local 
time [after Mead and Fairfield, 1975]. 


A simplified version of the eircurt 
shown in Figure 62. Cu = Cy + C- 


This circuit determines the frequency 
of the LC-oscillations thought to 
represent the Pc 5 micropulsations. 
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FIGURE PAGE 


66 Latitude profiles of magnetic perturbations 160 
in geomagnetic coordinates. Each successive 
profile is-separated by ~ 2 minutes. The 
profddespan. the, left columns (i,ec4 profiles 
A, C and E) show that the AY step occurred 

at low latitude while the profiles in the 
right column (i.e., profiles B, D and F) 
Show that the AY step occurred at higher 
latitudes. The borders of the westward 
electrojet are indicated by arrow heads at 
the bottom of each frame. Note the motion 
of the equatorward border. 


67 A theoretical latitude profile along a 162 

| meridian 4° to the west of the central 
meridian of the ionospheric portion of the 
model current system. The model three- 
dimensional current system involves both 
the westward Hall current, southward 
Pedersen current and Birkeland sheet 
Currents. The current system has a 
latitudinal extent of 4° and longitudinal 
extencaotse | 6428 Aycurrent of 10° A was 
used in the calculation. 


68 A differential profile for one quarter 164 
cycle of the oscillation in borders and 
intensity of the current system shown in 
Figure 67. 


69 Two differential profiles one quarter of a 165 
cycle apart. The data were bandpass fil- 
tered (4-100 mHz). The perturbation 
pattern was associated with micropulsations 
Athy noamhzZe 


70 Magnetograms for the Day 347, 1971 event. V2 

7a 1-20 mHz bandpass filtered magnetograms 173 
Foretnendayss475 LOZ Leavenk 

ic Magnetograms for the Day 253, 19/74 event. 174 

13 1-20 mHz bandpass filtered magnetograms 175 


forathe : Day .253,.1974<event. 
74 Magnetograms for the Day 254, 1974 event. [als 
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FIGURE PAGE 


75 1-20 mHz bandpass filtered magnetograms V77 
for the Day 254, 1974 event. 

76 Magnetograms for the Day 255, 1974 event. EZS 

77 1-20 mHz bandpass filtered magnetograms 179 


for: the Day 255, -1974 event. 


78 Raw power spectra for the three components 182 
CHenDeucyetor the Dayns47e41977 (1745- 
1845 UT) event. The data were filtered 
with a 1-20 mHz bandpass filter. 


79 The locations of the convection westward 184 
electrojet plotted as a function of 
invariant latitude and Universal Time 
for the Day 347, 1971 event. 


80 Same as Figure 78 except for the Day 253, io5 
1974 (1330-1430 UT) event. 


8] Same as Figure 78 except for the Day 254, 186 
1974 (0945-1045 UT) event. 


82 Same as Figure 78 except for the Day 255, 188 
1974 (1000-1100 UT) event. 


83 Latitude profiles of power and relative 189 
phase along with the polarization ellipses 
for the Day 347, 1971 (1745-1845 UT) event 
attdl e2umAzecothee locationscofuthers tations 
are indicated at the top of the figure. 
The ellipses that are shaded indicate CW 
polarization and those that are not shaded 
indicate CC polarization. 


84 Longitude profile showing variations of 19) 

power, relative phase and polarization. 

The longitude is measured in the geomag- 

netic coordinate system. MCMU is % 3° 
equatorward of the three other stations 

(which lie approximately along the same 

latitude line) and is due south of URAN. 
lineepdot. 1s for the Day 253, 19/740(1330- 

1430 UT) event at 12.5 miz. 
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FIGURE PAGE 


85 Same as Figure 84 except for the Day 254, 192 
1974 (0945-1045 UT) event at 11.7 mHz. 

86 . Same as Figure 84 except for the Day 255, Ls: 
1974 (1000-1100 UT) event at 11.7 mHz. 

Al Block diagram of the basic magnetic recor- Ziae 
ding system at each station. 

x2 Gain curve of the magnetometer. 224 

A3 Theoretical frequency response curve for (e(4|s) 

‘ the magnetometer. 

A4 Block diagram of the digital recording Ciel, 
system. 

Bl The spectral window for smoothing the (2 8)/i 
periodograms. 

B2 Variation of the confidence interval and 238 
bandwidth. 

Gl A latitude profile for an E-W current system 241 


with 20° E-W extent [after Kisabeth, 1972]. 

C2 Latitude profiles associated with current 243 
Systems having various latitudinal current 
distributions [fafter Kisabeth, 1972]. 


C3 A latitude profile for a N-S current system 245 
with 4° E-W extent [after Kisabeth, 1972]. 


C4 A latitude profile for a N-S current system 246 
with, 4° E-W extent [after Kisabeth, 1972]. 
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CHAPTER I 
INTRODUCTION 


1.1 Definition and Classification of Micropulsations 
| The earth's magnetic field is not constant but 
fluctuates with time with periods ranging, from small frac- 
tions of a second to more than 30 millions of years. While 
the longer period changes are of internal origin, variations 
with periods less than a few days are of external origin. 
That part of the geomagnetic spectrum which covers a fre- 
quency range of 1.7 to 5000 mHz with the corresponding period 
range of 0.2 seconds to 10 minutes is generally referred to as 
geomagnetic micropulsations. Their amplitudes can vary from 
a fraction of a gamma to hundreds of gammas (1 gamma = 
1ont = 10° gauss). This thesis deals with micropulsations 
in the 1 to 16 mHz frequency range. 

The history of micropulsations stretches back over 
more than a century. In 1861, Stewart, while studying a 
great magnetic storm recorded at Kew Observatory, noted: "The 
exertion of the disturbing force was of a throbbing or pulsa- 
tory character. The interval of time between two of these 
Rint TeEeaone may be said to have varied from half a 
minute, or the smallest observable portion of time, up to 
four or five minutes." [Stewart, 1861]. Since then, a great 


diversity of papers has appeared in the literature dealing 
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with the topic of micropulsations. The largest increase in 
the output of papers came about in the late 1950's after the 
I.G.Y. More recent advances in the study of pulsations 

have stemmed from the improvement of instrumentation, viz. 
the development of induction magnetometers, establishment of 
magnetic Pape yecording Systems, etc. 

Studies of the characteristics of micropulsations 
have led to the separation of this activity into various 
groups and many authors have devised their own classifica- 
tion systems. In order to avoid confusion and to standardize 
the nomenclature used, a classification system was approved 
at the thirteenth General Assembly of the IUGG in Berkeley, 
Gahifornia iidacobs et-al., 1964]. In this’ classification 
System, micropulsations were divided into two main groupings: 
those of regular and continuous appearance called Pc, and 
those with an irregular form called Pi. The Pc group was 
further divided into five subgroups and the Pi group into 
two subgroups. The classification is shown in Table 1. It 
should be emphasized that the properties of different classes 
may overlap and that there may be different source mechanisms 
producing pulsations within the same class. Until a better 
classification system is devised based on the correlative and 
genetic properties of micropulsations, the above classifica- 
tion based on morphological properties will continue to be 
useful. This thesis is concerned with the regular and con- 


tinuous pulsations in the spectral range 1 - 16 mHz in the 
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morning sector. Thus, the pulsations discussed in this 


thesis belong to the groups Pc 4 and Pc 5. 


1.2 The Magnetosphere 

Since micropulsations are of external origin, a 
brief description of the earth's magnetic field and plasma 
environment is warranted here. 

The earth is enclosed in a cavity, called the 
magnetosphere, which is carved out of the solar wind by the 
earth's magnetic field. The basic configuration of the 
magnetosphere is shown in Figure 1. 

The idea that a stream of particles flows outward 
from the sun into interplanetary space originated as long as 
300 years ago with the work of Halley. Near the beginning of 
this century, Stoérmer [1906] and Birkeland [1908, 1913] intro- 
duced the concept that charged particles of solar origin 
might be responsible for the production of geomagnetic storms 
and auroras. Chapman and Ferraro [1931] described the solar 
plasma from the hydromagnetic point of view in the early 
thirties. Their model suggested that discrete clouds of 
plasma were emitted from the sun during solar flares. Guided 
by the observation that the type I tails emitted by comets 
(tails of ionized material which scatter sunlight) always 
point away from the sun regardless of the position of the 
comet on its trajectory, Biermann [1951] developed a more 


realistic picture which featured continual flow of plasma away 


= 


tae. ni beeeus2ib znoltsetug ‘oath eeu aba 
) 18 3% pag 4 quay sit? ot 7 


_ proneacdoveat 
& intptys ‘engages ¥O St 2nofde2tuqoxo tm bing sot 


pea 

-< 

ao 
7] 


ef brs bIstt si deqpsm ot ites ‘eng 0 motte ' i 
sser Reo au 


BAe bol feo Wives ra ti bezelfons zt ‘Ao, itt 


ant vd boarw Yaloz add to yo bavitsD Bt Hotline 


* 


gftt Ho 10 i¢ sue Finds ‘oteed SAT blot on E 


mal, i 


| ovup Ha ni none) ah 


DiBWILO Zwolt 2afott veg. 7.0 magnet 2 f S683. seit, 9 


hows 


fof 26 bstsnatptio 9D 6GR Avot oasiiqres at, tere nt 


70 prifnfged aid 165i Nob Ustt 70 sow Sad cee ‘4 | 
-ovsni TEST | eae Fi bas totais bin rade ry Sar 2g rer 
| ntpty+o Veter to eslotiveq begisno. ter saps $ 


anitbaspa | 2Ftengsmosy to nots asibo-a: out 43 ‘oT ay 


7 t vou 


vsloe ont bsdisa2eb reed] ati. bn. gna ; atom 


nt a. 
| yfyss Slt ne woty. +s Intoy > FF aivpainoy ryt ant ort ue 
ee nrag 
+0 zbuoto. stensath. Sarid besesopue fabom whant 2 “oe Pd 
? 5 if i ie, : 
ms 7 


Tq. 


bebing eave? Bares aataab! me ont mont batttas anew: 
2Fsno3 yd. beitine et isa I Sgvs sad tend nota svasede 
eyew!s’ (sip Fniie 4eijoe' 4 salt dw ietyatem basinat oe i 
sad yo nots tzog. aad. t6 zesTbnseot fiv2 sat pont win tig | 
10m 6 ‘bogotavab: [ree sy: nem th “\sotoeters att a0 3 Si | 


ee % 4 a 
us ‘amas tq Fomor feyotsaes” boruteat aorde- santo tq ot 


M 


a is ym 


| aunbe4 


HS 
MOG 


PSSA OSECCCY.O 
SOOKE OOOO e OO ae 
etna eeronetataneeaoneracenenet O55 revo 


° are wR 
e ges el 
a SA. 
oe 4 > ee 
tas ; ~ Par a PH a 
Sie ‘SaeRA Igy ey em: : 
S és VAs; M4 EA 5 
< 2 ig A, om 
~ > fi} ald, fi: of $ 
~ SS a Pfs KGS a . 
e % 
e f § 
vi S = i ate & 
e Le Meee 
4 @ of.) 8 

“I . at .*, ei eo © 

<a ae et ee ty o@ ' ° Z 
0” of =A lo 
i ae. ay” 
puss ok : 
og Soh 
er ere nie 
° 


° q °P 
Pop.-to doy e 


HONOUL 


SSITIILYVd 
~ QaddvVeL 


(dSND) 
~14]319 


~~ 


f sEF4 


6 


from the sun (which is now known as the solar wind). After 
Biermann's contribution, Chapman [1959] in the late fifties, 
developed a model for the solar atmosphere in static. 
equilibrium and predicted that a finite pressure must be 
exerted inward on the solar atmosphere at great distances. 
Parker [1960] pointed out this error of finite pressure at 
infinity in Chapman's theory. In a series of papers beginning 
in 1957 on the investigation of the expansion of the solar 
corona into interplanetary space, Parker developed a general 
theory of steady, radial and spherically symmetric flow of 
plasma from the sun. [See Parker, 1963]. 

The investigation of interplanetary Space by arti- 
ficial satellites and space probes has confirmed the model of 
the solar wind as developed by Parker. It is now established 
that the solar -wind,is a completely ionized medium composed 
primarily of hydrogen (protons and electrons) with about 5% 
of ionized He streaming continually from the sun. The solar 
wind carries with it the coronal magnetic field from the sun 
which is termed the interplanetary magnetic field (IMF). 
Approximately 80% of the IMF is confined to the ecliptic plane 
with a small component perpendicular to the ecliptic plane. 
This normal component of the IMF seems to play a major role 
in the initiation of magnetospheric storm and substorm activi- 
ty [Dungey, 1953, 1961; Fairfield and Cahill, 1966; Rostoker 
and Filthammer, 1967]. Substorm activity tends to be enhanced 


less than an hour after the normal component of the IMF has 
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turned southward (i.e., pointing below the ecliptic plane) 
[Arnoldy, 1971; Rostoker et atS5 97258 Meng! et alts4 9739¢e 
Typical values for the solar wind parameters at the 
orbit of the earth are: bulk velocity ~ 300-500 km/sec, num- 
ber density 5-10/cm>, , IMF ~ 5-10y (the field strength 
4 


increases with high solar activity), temperature ~ 5x10" °K 


(It can vary between 10° °K and 10° 


°K. The temperature tends 
to be higher when the sun is active.), Alfvén velocity 
~ 50-100 km/sec and velocity of sound in the solar wind 
~ 100-200 km/sec. 

Since the solar wind is both supersonic and super- 
Alfvénic, a standing magnetohydrodynamic shock wave is set up 
approximately 14 Re in front of the earth. Downstream from 
the standing shock, the bulk flow becomes subsonic (the Mach 
number being ~ 0.2) and the solar wind plasma is thermalized. 
This region of shocked plasma is called the magnetosheath. 
The position at which the dynamic pressure of the shocked 
solar wind plasma is balanced by the outward pressure of 
earth's magnetic field marks the boundary of the magnetosphere, 
called the magnetopause (which is about 100 km thick). On 
the dayside, the magnetospheric magnetic field lines are more 
or less dipolar and have a magnitude of ~ 50y just inside the 
magnetopause. However, under the sweeping action of the solar 
wind, the magnetic field lines are greatly distended on the 
night side forming a magnetotail which reaches beyond the 


orbit of the moon at ~ 60 Re- The center of: the tait «1s 
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8 
dominated by a region of energetic plasma (number density 
% 0.1-30/cm?) known as the plasma sheet. The plasma sheet 
thickness increases from about 5 Re near* the center of the 
tail to about 10 Re near the flanks. Its width is. about 20 
Re. The center of the plasma sheet marks a region of magne- 
tic field reversal known as the neutral sheet (which is about 
600 km thick). -On either side of the plasma sheet are re- 
gions of low plasma density (< 0.01/cm?) and enhanced magnetic 
field termed the tail lobes. Typical magnetic field strengths 
in the near earth plasma sheet are ~ 10y while the tail lobe 
at similar distances from the earth has field strength of 
15-20y. 

The magnetosphere is populated by 'cold' plasma, 
'hot' plasma and energetic electrons and protons. The ‘cold' 
or thermal plasma refers to particles with energies in the 
range of a few eV to a few tens of eV. They form a region of 
dense plasma surrounding the earth called the plasmasphere. 
This volume of cold plasma has a toroidal shape with cross 
section resembling the dipole field. The plasmasphere has a 
field aligned boundary called the plasmapause. The plasma- 
pause separates the region of high thermal plasma density, the 


3 in the equatorial plane) 


plasmasphere (~ 100 electrons/cm 
from a region of low plasma density, the plasma trough 
(v ] electron/cm* in the equatorial plane). The mean equator- 


ial radius of the plasmapause is typically about 4 Re» but may 
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vary from about 7 Re during extreme quiet periods to 2-3 Re 
during disturbed times. Outside the plasmapause, the thermal 
plasma density decreases approximately as Ri Inside the 
plasmapause, the density distribution along field lines is 
more or less in agreement with a diffusive equilibrium model 


4 model. 


which gives a much less rapid decrease than the R- 

The ‘hot' plasma refers to particles whose energies 
range from several tens of eV to several tens of keV. This 
plasma population, though less dense than the thermal plasma, 
carries the largest share of the particle energy density. The 
energetic particles are particles with energies that range 
from a few hundred keV to a few MeV (for electrons) and a few 
hundred MeV (for protons). They are the least numerous group 
of particles and constitute the so-called trapped Van Allen 
radiation belts. | | 

It is now well known that energy from the solar wind 
penetrates to the interior of the magnetosphere and is trans- 
ferred from one place to another within the magnetosphere 
through the convective motion of the flux tubes (and their 


> 
frozen-in plasma) Wis = si as well as the gradient drift 
B 


& We 2 : 
[Vo cra; Boxe VBA. and curvature dri fc 
eB 


| MKS units are used. E is the electric field and B the 


magnetic field. 
ae = 5 mv is the K.E. of the particle associated with the 
gyration around field lines. e is the charge of the particle. 
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Recent studies by Rosenbauer et al. [1975] suggest 


(B- °V By]? of the energetic particles. 


that solar wind plasma gains entry to the magnetosphere 
through the polar cleft and becomes part of the plasma mantle 
streaming downtail away from the earth. It is thought that 
the plasma mantle particles eventually populate the plasma 
Sheet and acquire a drift velocity and field-aligned stream- 
ing velocity toward the earth [Rostoker, 1976]. 

Insofar as the motion of the plasma sheet particles 
close to the earth is concerned, curvature and aradi ene drift 
become significant causing electrons to drift eastward and 
protons westward. The drifting particles have a component of 
their population which precipitates into the upper atmosphere 
both acing and ionizing the constituent atmospheric parti- 
cles. It is now understood that the plasma sheet of the 
magnetotail maps into a region of enhanced auroral luminosity 
and electrical conductivity called the auroral oval [Frank, 
1971: Lassen.) 19744eCbicure 2). An artustesaconception of 
the auroral oval is shown in Figure 3. 

Rostoker and Bostrém [1976] have suggested that the 
electric field driving the ionospheric current originates 
in the magnetosphere as shown in Figure 4. Since the plasma 
in the magnetotail is collisionless, one may consider the 


conductivity along the magnetic field lines to be nearly 
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infinite, in which case the magnetic field lines may be 
considered as equipotentials. Thus any electric field normal 
to the magnetic field in the magnetosphere will map along 
magnetic lines of force into the ionosphere. The magneto- 
Spheric electric field when mapped into the ionosphere is 
Shown in Figure 5. As mentioned earlier, the auroral oval is 
a region of enhanced conductivity due to energetic particle 
precipitation. Thus, the existence of electric fields in a 
region of high conductivity results in electric current flow. 
In Figure 5, it can be seen that the Hall currents 
(current flow whose direction is perpendicular. to both the 
electric field and magnetic field) are driven in the auroral 
oval such that the flow is eastward in the evening sector (the 
eastward electrojet) and westward in the morning sector (the 
westward electrojet). Based on the electric field data of 
Heppner [1972] and Mozer and Lucht [1974], it is considered 
that the above mentioned electrojets are Hall currents flowing 
in the E-region of the ionosphere (~ 100 km above the surface 
Ofethe earth). in addition. to. the: Hall currents, direct 
(Pedersen)* currents tlhow in: the directronm of the sn=seelectric 
fieldy) ihe.conductivity discontinuities) at the edges or the 
auroral oval cause the ionospheric current to diverge and flow 
along the magnetic lines of force. These field-aligned 
currents (the Birkeland currents) flow into and out of the 


tail region and couple the magnetosphere to the ionosphere. 
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They have been observed using magnetometers aboard polar 
orbiting satellites [Zmuda and Armstrong, Vee and are 
thought to be driven by MHD generation processes in the 
magnetotail [Rostoker and Bostrém, 1976]. The flow pattern 
of these currents is already shown in Figure 2. Finally, we 
note that there tends to be net downward current flow in the 
morning sector and net upward flow in the evening sector 
[Yasuhara et al., 1975]. As we shall demonstrate later in 
this thesis, Pc 5 micropulsations in the morning sector which 
have obatt amplitudes in the auroral oval are generated by a 
combination of spatial and temporal oscillations of the 


current systems described above. 


i337 Ine Morphology of Pc jis andiw Pct4' 

(aseP’c 5 

This class of pulsations with the longest period 

has amplitudes of a few gammas to several hundred gammas. The 
oscillations are rather regular, with a wave train consisting 
of a series of damped type variations lasting about an hour or 
more. An example of Pc 5 recorded at high latitudes is shown 
in Figure 6. The duration and amplitude of the wave trains in 
the morning hours are generally longer and greater respectively 
than those in the afternoon. The oscillations are most sinu- 
soidal in the morning and early afternoon [Jacobs and Wright, 
1965]. The Z-component (the vertical component) of such 


pulsations can reach appreciable amplitude in the auroral 
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18 
zone [Eleman, 1966] or at the latitude of maximum intensity 
of the horizontal components [Samson, VI7212" Argiven Pc=5 
event has the same frequency at all latitudes [Ellis, 1960; 
Obertz and Raspopov, 1968; Samson and Rostoker, 1972]. There 
are only a few exceptional cases in which the H-component 
varies with latitude [Siebert, 1964; Voelker, 1968; Rostoker 
and Samson, 1972]. Oscillations similar in amplitude and 
period often appear simultaneously at conjugate stations 
[Nagata et al., 1963; Jacobs and Wright, 1965]. 

It is now well established that the occurrence 
frequencies and amplitudes of Pc 5 maximize in the auroral 
zone (65° - 70° geomagnetic latitude) [Jacobs and Sinno, 
1960; Kato and Saito, 1962; 01', 1963; Kokubun and Nagata, 
1965; Obertz and Raspopov, 1968; Hirasawa, 1970; Gupta, 1973, 
1974]. Samson [1972] found that the intensity maximum of 
Pc 5 follows the statistical auroral oval as illustrated in 
Figure 7. 

Gupta [1973] carried out a statistical study of the 
occurrence frequencies of Pc 5 using seven Canadian stations 
covering a latitudinal range (83.1°N - 54.3°N) and a longi- 
tudinal range (28/.7°E -— 347.2°E). He: round that the diarnal 
variation in occurrence frequencies of Pc 5 shows a peak in 
the morning and a peak in the evening hours for all stations. 
The evening peak is more dominant in the polar zone. In the 


auroral zone, the morning peak is more dominant in agreement 
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with the results of Saito [1964] and Hirasawa [1970] who 


noticed that the occurrence frequency peak in the morning is 
Sharply defined and centered near 0400-1000 LMT. Gupta also 
noted a midnight peak in the auroral zone. In the sub-auroral 
zone, the evening peak is more pronounced. Saito [1964] noted 
that Pc 5 are observed with equal probability over most of 
the day in mid latitudes. The basic diurnal variation pat- 
tern of the occurrence frequency of Pc 5 is unchanged for 
different seasons and different levels of magnetic activity 
[Gupta, 1973]. While simultaneous occurrences of Pc 5 rarely 
were noted at the seven Canadian stations (about 2% of the 
events), Gupta found simultaneous occurrence in the auroral 
zone stations covering an area of 10° latitude and 20° in 
longitude in more than half of the events he studied. Obertz 
and Raspopov [1968] observed that the characteristic longi- 
tudinal extent of Pc 5 is 60°. The following results .of 
Gupta [1973] should be noted in passing: 27 day recurrence 
tendency is weak; correlation between occurrence and Kp is 
dood,up toyKp =.5° and Kp -=0373 correblationebetween 
occurrence and daily and monthly mean relative sunspot numbers 
is poor; there is greater occurrence of Pc 5 in summer than 
in winter. 7 

The diurnal variation in amplitude of Pc 5 shows a 
peak in the early morning and another in the evening [Jacobs 
and Sinno, 1960; Kato and Sato, 1962; Kokubun and Nagata, 


1965] with the morning peak dominating the evening one. 
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Gupta [1974] noticed that, in the auroral zone, the diurnal 


amplitude maximizes in the morning and this peak tends to 
move northward reaching the polar station of Resolute Bay 
around noon. Saito [1964] also observed a single maximum in 
the diurnal variation of amplitude at the high latitude 
stations of Fort Churchill and Point Barrow. Only occasion- 
ally did Gupta find a maximum of amplitude during prenoon 

and post noon hours for the seven Canadian stations. He also 
reported a midnight peak for all the auroral zone stations and 
that the diurnal variation pattern of Pc 5 amplitude showed 
no clear relationship with seasons or with magnetic acura. 
The amplitudes of Pc 5 do not show 27 day recurrence tendency 
but seem to have a linear relationship with Kp up to a level 
of 6- 

The polarization of Pc 5 in the horizontal plane 
exhibits: both latitudinal and diurnal’ variations. Kaneda et 
al. [1964] and Obertz and Raspopov [1968] showed that Pc 5 
exhibit pronounced latitudinal change in the sense of polari- 
zation. They reported that Pc 5 at Barrow (68.6°) had cw 
polarization (clockwise polarization looking along the direc- 
tion of the magnetic field) and Pc 5 at College (64.7°) and 
Sitka (60°) had cc (counter-clockwise) polarization when the 
center of activity was between College and Barrow. An 
appreciable latitudinal phase shift in the H-component 


accompanied the polarization reversal. Samson et al. [1971] 
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also found that Pc 5 observed at stations south of Fort Smith 
(67.3°) had a sense of polarization opposite to that at 
Cambridge Bay (76.8°). In addition, Samson et al. [1971] 
noted that the line of intensity maximum of Pc 5 coincided 
with the demarcation line separating the regions of opposite 
polarization. 

Earlier, before the discovery of the demarcation 
line, studies of the diurnal variation of Pc 5 polarization 
in the northern auroral zone region by Nagata et al. [1963], 
Kato and Utsumi [1964] and Kokubun and Nagata [1965] indicated 
that the sense of polarization in the horizontal plane was 
predominantly cc in the morning and cw in the evening. Samson 
et al. [1971] also observed the above polarization patterns 
equatorward of the demarcation line. A summary of both the 
latitudinal and Adena variations of the H-D polarization of 
Per 5 ise shown in Figure 8. 

The polarization characteristicswor i mroad? band oft 
micropulsations have been investigated statistically by Rankin 
and Kurtz [1970]. They have analysed the magnetic data 
recorded at 14 stations distributed in central and southwest 
Alberta from 49° 51'N geographic to 54° 26'N and 112° 17'W 
geographic to 117° 50'W. The average senses of polarization 
in the horizontal plane of micropulsations in the spectral 
range 3.3 - 100 mHz are shown in Figure 9. This spectral 
range covers the higher frequency bands of Pc 5, and the 


complete Pc 4 and Pc 3 bands. It can be seen from Figure 9 
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25 
that the diurnal polarization pattern of high latitude Pc 4 
and Pc 5 pulsations are similar, while the sense of polariza- 
tion of Pc 3's (22.2 - 100 mHz corresponding a period range 
of 10-45 seconds) is opposite to that of Pc 4, 5's with cw 
polarization in the morning and cc polarization in the after- 
noon. 

The cc polarization in the morning and cw polariza- 
tion in the evening are valid not only for the statistical 
case but also for simultaneous observations on both morning 
and evening sides of the earth [Troitskaya, 1967]. This 
characteristic in polarization has led Nagata et al. [1963] 
and Kokubun and Nagata [1965] to suggest that Pc 5 pulsations 
are due to transverse hydromagnetic waves on magnetic field 
lines near the magnetopause excited by a hydromagnetic inter- 
action of the solar wind and the geomagnetic field. Atkinson 
and Watanabe [1966] subsequently suggested that the transverse 
hydromagnetic waves are a result of a surface wave motion at 
the magnetopause and the surface wave is generated through 
the development of Kelvin-Helmholtz instability as first 
suggested by Dungey [1954]. Based on this distinct diurnal 
variation pattern of polarization in H-D plane, Samson et al. 
[1971] also tended to favour the idea that Pc 5 pulsations are 
generated through the development of Kelvin-Helmholtz in- 
stabilities at the magnetopause. A more detailed discussion 


on the theory of micropulsations will be given in the next 


chapter. 
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The diurnal variation of the sense of polarization 

is reversed in the southern auroral zone compared to the 
northern auroral zone [Nagata et al., 1963]. This is to be 
ane Since the senses of polarization are conserved along 
field lines [Nagata et al., 1963; Jacobs and Wright, 1965]. 
In mid latitudes the diurnal polarization characteristics 
are more complicated with four sectors alternating from cw 
BOMCC  oaito.,7 1964 |. 

, With regard to the polarization ellipses in the 
H-D plane, Samson [1972] observed that the morning sector 
has consistently elliptical cc polarization and that the 
afternoon sector has a mixture of senses of polarization and 
a predominance of almost linear polarization in H. The 
ellipses at the high latitude stations (58.4° - 76.8°) are 
predominantly polarized in the H direction. 

There have been only a few studies on the polariza- 
tion characteristics in the vertical planes for a few events 
[Paulson et al., 1965; Wilson, 1966; Annextad and Wilson, 
1968]. Samson [1972] carried out a statistical study on the 
polarization characteristics in the H=ZvandQ0=7" planessn Nos t 
of the stations show cw polarization in the H-Z plane over the 
entire day. The ellipses are predominantly polarized in the 
H direction. From the analysis of the polarization character- 
istics in the H-Z plane, Samson concluded that the perturba- 
tion magnetic field of Pc 5 is definitely not polarized trans- 


verse to the main geomagnetic field. The senses of polariza- 
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tion in the D-Z plane at Cambridge Bay (76.8°) and Fort Smith 


(67.3°) are cc in the morning and slightly cw in the after- 
noon and evening with the reverse being true at the southern 
Stations. The major axes at most stations seem to favour the 
D direction in the interval (5ic5 29 cl ain eehem-Z olanes, At 
other times, the polarizations are more confused. 

This thesis is concerned with the Pc 5 pulsations 
in the morning sector which have '‘cleaner' wave forms and 
‘clearer’ polarization characteristics than pulsations 
occurring at other times of day. The Pc 5 pulsations which 
occur during the main shag of geomagnetic storms and are 
confined to the afternoon sector [Barfield et al., 1972; 
Lanzerotti et al., 1975] and the midnight Pc 5 noted by 
Gupta [1973, 1974] will not be dealt with in this thesis 
since these pulsations, though in the Pc 5 range, are 
probably due to different origins. 

Gb)'e Shc 24! 

Pc 4 is described as a regular and continuous 
pulsation with period between 45 and 150 sec. This existing 
Glass huication.fails.to,describe,. the unigquerjeatumes of sche 
kind of pulsations reported by Rolf [1931], Harang [1932, 
1936..1939,.1940-1944]. Sucksdorff. 11939], Veldkamp,[ 1960] 
and Annexstad and Wilson [1968]. They referred to that 
particular pulsation within a narrow band in the Pc 4 range 
as giant pulsations or Pg's. However these are not giant 


pulsations in the literal sense because their maximum ampli - 
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tudes only reach 30 gammas or so which is comparatively much 
less than those of Pc 5 which reach hundreds of gammas. We 

Winrenter to this unique pulsation as “Pc 4° to be distin- 

guished from other Pc 4's such as those associated with the 

plasmapause studied by Lanzerotti et al. [1974]. 

Atypical example of'Pe 4° yecorded @tour sta- 
tions is shown tn Figure 10. It can be seen that the wave 
form is highly sinusoidal and well modulated. The average 
period of 'Pc 4' is 87 sec and its duration varies between 10 
minutes and several hours [Sucksdorff, 1939]. 

Thesoccurrence of ‘Pcl 4°:ti's? rare. Eb Rolf [1931] 
reported 28 occurrences from almost 9 years of five from the 
auroral zone. station of Abisko. Sucksdorff [1939] observed 
150 events over 25 years (1914-1938) of data at Sodankyla 
with an average of 6 ‘Pc 4' events per year. Korobkova et al. 
[1959] found not one single case in the 20 years of data at 
TikhayasBay (71.1°) and none were obServed at the Soviet 
polar Stations during 1.G.Y¥. Rolf [1931], Harang [1932, 
1940-1944) and Sueksdorff)[1939), noted thatuerc’ 4 occurs 
most frequently during early morning hours (a-7' LT). around 
equinoxes and in years of minimum solar activity. The ‘Pc 4! 
pulsations also tend to recur for successive days at intervals 
of about 24 hours [Sucksdorff, 1939]. 

The 'Pc 4' pulsations are extremely localized in 


latitude. This feature is well illustrated by an event on 
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July 17, 1958 studied by Veldkamp [1960] using a large number 


of stations in Europe. The latitudinal distribution of ampli- 
tude of that particular event maximized at ~ 57°N geomagnetic 
and the amplitude dropped to the noise level about 5° north 
and 5° south of the latitude of maximum intensity. The 
pulsations occurred in a strip confined within 50° to 60° 

(The overall spatial distribution of the pulsations was 
confined to a small area about 1000 km by 1000 km.). The 

times of beginning and ending of the pulsations were differ- 
ent for different stations and almost simultaneous occurrence 
was restricted to adjacent stations which were less than 500 
km apart. Most stations registered cc polarization in the 
horizontal plane. Veldkamp contended that the observed period 
of 105 sec fits well with eigen oscillations of period ~ 100 
sec predicted by Obayashi and Jacobs [1958] for field lines 
emanating from 57° geomagnetic latitude. 

Annextad and Wilson [1968] analysed 20 'Pc 4' events 
recorded at College and its conjugate, Macquarie Island during 
years of minimum solar activity. None of the 20 events was 
Seen at polar and mid latitude stations. /5, of the events at 
College were polarized in cc sense on horizontal plane while 
67% of the events at Macquarie Island were polarized in cw 
sense. Therefore, continuity of the sense of polarization was 
preserved along the magnetic field lines. 70% of the ellipses 


in the horizontal plane had their major axis aligned within 
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Beto TUCO POR=Wh) This medns Sthat 'Pc 4° has a larger D-component 


than the H-component. The polarization ellipses in the H-D 
plane for the two stations were mirror images of each other 
across an E-W line. Annextad and Wilson contended that the. 
‘Pe 4' is a result of the n = 2 (first even tode) mode of a 
transverse hydromagnetic wave propagating along field lines 
to conjugate points. 

It should be noted that Harang [1939] noticed 
Simultaneous occurrence of pulsations on magnetic records and 
in the scattering region of radio waves at heights of 650-800 
km. He used a powerful transmitter at a frequency greater 
than the penetration frequency Foo (i.e., the minimum fre- 
quency at which the O wave can penetrate the F layer) and 
obtained pulsating echoes above the F layer. These pulsating 
echoes coincided both in duration and period with 'Pc 4! 
pulsations. 

We were fortunate to observe a few occurrences of 
Gheselrare Pc A" pulsating in our'idata andy tney-are docu - 


mented in this thesis. 


A more complete review on micropulsations can be 
found in the following reviews: Jacobs and Westphal [1964], 
Hultqvist [1966], Troitskaya [1967], Campbell [1967], Saito 


[1969], Jacobs [1970] and Orr [1973]. 
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1.4 Purpose of This Thesis 


A total of nine magnetic stations were in operation 
in 1971 and early 1972 in Western Canada, arranged approxi- 
mately along a corrected geomagnetic meridian, and providing 
very good spatial coverage for high latitude micropulsation 
events. In particular, the close spacing of stations in the 
auroral zone provides reliable information on the motion of 
the convection auroral electrojet. 

It is the object of this thesis to present the 
relationship between the Pc 5 micropulsations and the convec- 
tion westward electrojet and to demonstrate that the genera- 
tion mechanism for Pc 5 is due to the oscillation of the 
three dimensional current system associated with the convec- 
tion westward electrojet. 

In addition, an analysis and interpretation of the 
rare 'Pc 4' micropulsations recorded during 1971 and 1974 
(when three of the stations were distributed along the same 
latitude line) will be presented, and we shall show that 
these events may well represent the pure resonance effects 
which have been proposed by several authors as a source for 


long period micropulsation activity. 
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CHAPTER II 
EXISTING THEORETICAL CONCEPTS 


Most of the attempts to explain pulsations in the 
earth's magnetic field are based on the concept of hydro- 
magnetic waves which were proposed by Alfvén in 1942 [Alfvén, 
1942]. This chapter begins with a brief review of the modes 
of propagation of hydromagnetic waves in a perfectly conduc- 
ting fluid permeated by a uniform magnetic field. This is 
then followed by a discussion of the modes of hydromagnetic 
oscillations of the medium in the presence of a non-uniform 
magnetic field such as a dipole field. The excitation 
mechanisms are discussed next, and the theories of the 
coupling of energy from an external source to a local field 
line are reviewed. The chapter is concluded by comments on 
the existing theoretical concepts as related to experimental 


observations. 


2.1 Hydromagnetic Waves 


A qualitative review of the subject of hydromagnetic 
waves is presented in this section. Much more comprehensive 
treatments can be found in any standard text on plasma physics 


[e.g., Alfvén and Falthammar, 1963]. 


The crudest approximation to plasma behaviour is to 


regard the plasma as a perfectly conducting fluid. In addi- 


33 


F oh isn SRE 


ardaanes. sna Teapioas are | 
, i 
ano pt snatapenog. atstaxs: oF ciecaamiiee eas to 120M 


“orb vr vo tqaone> ant nd bs asd aye, bia tt 2 ntvancam 
BUTT fAT SRO mi AavtTA yA pazonoyy s198 azine 2aveM.2 


=~ 


esbom snd to waiver tehic 5 At rm aniged retesna. ett ; 
-wbnos ¥idostvoq a ff esvew of Tone emetoNT To net st 
SE erat PIS ES abbenpem MUST rei iB xe bea nannay te 
DHIHKMO VY to e2ebon IAF to Oi aauoeit 5 ba vowtt 
NOT ihu-~don s Toa.sansestd on) ei. lt ah aft. Fo 
mortstions sn] pists oT oni é' eb dobe vay 9 

gid To 25ify00%d odi bas! fxn bezeaoeih: mp at 

- bist reef 94 g2'inde fat s3 x9 a6 iO: APIS Be 


to 2twogwos Yd bebuisnoa 2F eigBio ont. onion we 


| Patnantrsaxe of batKiay 26 atqsonoe, Taatapnoan 


aval atte 
SisanesmowbyA Fo teabdwe’ atid) to watvey. reves: a 
ovtenanerqmeo syom out. Snottsee erdd at Mere 
—— Bhizetq mo 3x97. by sbnsd 2 YAS af bao? 34 nga. 
feaer +t seis ct 183 bing “wy en 
eo 2t IwWoTvetsa | amasty ae hota satxorqgs taeburs oat : 
-tbbs al .btuft gattoubnos yisostieq-s 26 omaghquede 


34 


tion to the hydrodynamic Properties, the fluid possesses 
equally important electromagnetic properties. In order that 
strong interaction between electromagnetic and hydrodynamic 


phenomena exist, the inequality 


] 
Le. [ke ee 
Ho Bo 


must be satisfied, where L is the linear dimension of the 
conducting medium of eh thy, QO; permeability u, and 
conductivity o in the presence of a magnetic field B 
[Lundquist, 1952]. Under normal laboratory conditions this 
criterion is not satisfied, but the situation is very 
different in the magnetosphere where the linear dimension 
is quite large. 

Let us consider a nonviscous and compressible 
fluid of infinite conductivity permeated by a uniform 
magnetic field BO: Let us assume that B is curl free 
CAR es so Vie x Ee = 0), the magnitudes of the perturbed quantities 
are small compared with those in the unperturbed state and 
the displacement current is negligible compared to the 
conduction current. Then, combining Maxwell's equations 
and the basic equations of hydrodynamics, one can deduce that 
three distinctly different waves can be excited in the hydro- 
magnetic medium. They are the Alfvén waves (also known as 
the shear Alfvén waves or anisotropic Alfvén wave when 


propagating along field lines), the slow magneto-sonic wave 
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(or slow magneto-acoustic wave) and the fast magneto-sonic 
wave (also known as fast magneto-acoustic waves and, when 
propagating across field lines, modified Alfvén wave; 
compressional Alfvén wave or isotropic Alfvén wave). The 
wave frequency must be low compared to the ion cyclotron 
Frequency otherwise the ions will not partake fully in the 
wave motion. If the angle between the propagation direction 
and the field B is 6, then the phase velocity of the Alfvén 
wave iS equal to Va cos 8 and the phase velocities of the 
magneto-sonic waves are given by the two real and non- 


negative roots of the quadratic equation for the phase 


velocities v2 
4 ade adil 2 (a AW 2 a 
which yields 
2 2 4 4 a2 2 We 
5) eu (VA + Ve 4Vc Vp COs. tc) 


‘ pr « 
where Va @ elaayte is the Alfvén velocity, Ve - fe 1s the 


Ho? 
velocity of an ordinary sound wave, P is the plasma pressure 
and 9p is the mass density. 

The phase velocity Va cos 8 of the Alfvén wave is 
equivalent to motion of the phase fronts in the direction B 
with speed V,. The Alfvén wave is thus a transverse wave. 
The slow magneto-sonic wave propagates with a speed which is 


smaller than the lesser of Vo and Va. The fast magneto-sonic 
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wave propagates with a speed which exceeds the greater of Ve 
and Va: The two magneto-sonic waves are neither longitudinal 
nor transverse except when propagating along field lines or 
perpendicular to the field lines. 

These waves are most conveniently discussed by 
means of polar plots of the velocity vectors of the three 
waves as shown in Figure 11. The vector OV drawn from the 
origin 0 to a point V on any one curve represents the velocity 
Of a wave and also the direction of the wave normal. Let us 
consider the special cases of propagation along field lines 
and across field lines. 

It is found that a pure Alfvén wave with speed Va 
and a pure sound wave (ion acoustic wave) with speed Ve 
propagate in the direction of the magnetic field (i.e., 6 = 0). 
A simple physical picture can be obtained for the pure Alfvén 
wave as follows. Because of the infinite conductivity assump- 
tion, one can picture the magnetic lines of force ‘materia- 
lmized’ into elastic strings to. which, the .conductmg fiuid 
particles have been 'glued'. Thus, the magnetic lines of 


force possess inertia (provided by the ions) and a tension 


© When the fluid of density po is disturbed from rest, the 
lines of force will perform transverse oscillations (i.e., 
the direction of particle motion is perpendicular to the 
direction of propagation parallel to the undisturbed magnetic 


field ie the magnetic perturbation b is perpendicular to Bae 
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The phase velocity o : [ens ion 
p ve hai of the wave is equal to density 


(eee): Vp = a MD in analogy to the definitions of the velocity 
HAP 
Oo 2 


of the transverse wave propagating along a stretched string. 
whe transverse nature of this mode implies field Tine: 
Muidance of energy. The Alfvén waves cause no density 
perturbation of the medium and would have been excited had 
the medium been incompressible. This concept of transverse 
oscillations of magnetic lines of force was used by Obayashi 
and Jacobs [1958] to compute the plasma density distribution 
of the ‘outer atmosphere' making use of the observed period 
of micropulsations. The stretched string analogy enabled 
them to postulate that the characteristic period of oscilla- 
Bing. lines of. force, is; given Dy 4. - 25 » the integration 
being over the length of the magnetic line of force. 

The sound wave that also propagates parallel to B 
involves no electromagnetic effects (in contrast to the 
Alfvén wave which is a purely hydromagnetic phenomenon) and 
is. a, longitudinal.wave. (i.e... particle motion isaparallel to 
the direction of propagation which is parallel to Ba) This 
mode can be thought of as a sound wave travelling in a pipe, 
the rigidity of the pipe corresponding to the stiffness of the 
magnetic field. 

There is only one mode that propagates perpendicular 
to the field lines (i.e., 9 = 90°). This is the fast magne- 


to-sonic wave moving across the field with a speed given by 


party fe badoa S482 pAOTS entdapaqoy pvew s2ieve , ‘3 
weet Spett ear ign! sbom Beg V9 Baniien sciaven 
obese On geueo pevew nav Tae aT 1BA8.19 cre 

ber betioxa 789d sved. bluew bing nis Poo aid +o wait a 
get avensid to 2g 53n00 eat oFdtezsvgqmagpt. fised ¢ 16a 

raz SYedo vd bo2u ao Bae TO. 2ente Asadeite 70 209 

woheudiiterb, ytteasb: leas a lay siuqnes oF paser] 


bo iraq bevyeedo os to gen pot dem STG eOnes vest? aa | 


. 


vk 


beidong Ypol sie pairyte hota sve git .2alrtee tut 
-sllio2d 30 botveg attziyedoe edd sid feds stnlessigie 
fic igvost tnt st. 9 La 1 vA nan ty ati es10t Yo ean 

| a70F FO, anti s+}50eam, Sat. 10: digest, + Nail 

of 98 [afieisq. zedepegong! oeté, ail svialy bia’ | q 
eit oF denrdnag nt) 2a26a7Fo! >ttampanonssate an, 
bie tnoinenonsti ot taagang bie Niptarg. Bei totter ci 

ot Tat eis at nottom efakaneg 2B * otew isa tg not 6-8 
oat al ae of harieneg, et ybaw Apttopeno-rg to. me r 
cant 6. tt pn travers ‘ow dbnvos is Et 10 tdpwoda pigs: a 
sith > evanttiae sift od! i PAoaeines aay tit 0 wet vey 
iad | = ey ay bier ‘ 
iit sonnet bom ano oe ai sa 
ecaad S26? ant Bt etd (Pe & ry yeaet) sant bts 
ee a esi scoean bhatt beh £26 as 


Ww 
7 : 
Pome. « an 
> 
er 8 r¢ 


: he i 


39 


ve + ve and causing compressions and rarefactions in the 
lines of force without changing their direction. This wave 
is essentially a longitudinal sound wave (i.e., direction of 
particle motion is parallel to the direction of propagation 
which is perpendicular to Bee B is parallel to B,) modified 


by the presence of the magnetic field. 


Cac POA Oidalitand To noi dale Osicidll ations 

Let us consider a hydromagnetic medium of infinite 
conductivity permeated by a nonuniform magnetic field B 
such as a dipole field. It is assumed that the medium is 
incompressible so that the acoustic mode waves in the plasma 
density are suppressed. Consider a thin shell surrounding 
the surface of revolution of a line of force. Then, upon 
being perturbed, the thin shell may experience two modes of 
oscillations of the medium; namely poloidal oscillations and 
toroidal oscillations. Poloidal oscillations can be 
visualized as compressions and expansions of the magnetic 
shells. The toroidal oscillations can be pictured as 
azimuthal twistings of the magnetic shells. Because of the 
nonuniformity of the magnetic field Bet extra terms involving 
the spatial derivatives of B are introduced into the 
governing equations. As a result, the poloidal and toroidal 
modes are coupled together. 

Dungey [1954] was the first to attempt to explain 


micropulsations in terms of poloidal) ands toroidalmoscit a= 
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40 
tions of the earth's magnetic field. Using Maxwell's 
equations and the fundamental equation of hydrodynamics, he 
was able to derive two coupled partial differential equations 
describing poloidal and toroidal oscillations. Since Dungey's 
work has had tremendous impact on the field of micropulsa- 
tions, it is worthwhile to see how the coupled partial 
differential equations are derived. 

Neglecting all external forces of non-electromag- 
netic origin such as those due to viscosity, gravity and 
pressure gradients and assuming the plasma velocity V to be 
small so that inertial terms can be neglected, the basic 


equation of hydrodynamics can be written in the form! 


ee (1) 


where 90 is the mass density of the plasma and dx Bis the 
mechanical force exerted by the magnetic field B on a volume 
element carrying the current density dieenaihe displacement 
current is negligible compared to the conduction current. 
Considering B as the sum of a constant dipole field 


B and a small disturbance b (160 >> (pale assuming Be to 


bepcurl. free (ise. , Vv x B = 0) and using Maxwell's equation 


vx B = Uads equation (1) becomes 


| MKS units are used. 
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Taking the vector product of the time derivative 
of equation (2) with BY and making use of Ohm's Law 
> 
Pee, =\ x B for a medium of infinite conductivity, equation 


(2) becomes 


Ee ge fg te ea ee (35) 
° > B 
Then, using the definition for Alfvén velocity Va = 172 
| (up) 
equation (3) becomes 
ehtai > 7 
Seat Sul Wein sco ee E) (4) 
ot 


Analogous wave equations for V and 6 may be similarly 
derived. | 

The above equation is the linearized hydromagnetic 
wave equation for an infinitely conducting fluid without 
plasma pressure gradients. It is the first stepping stone 
for many authors who have attempted to develop a theory for 
micropulsations [eg. Westphal and Jacobs, 1962; Orr and 
Mathew, 19713; among others]. 

Expressing equation (3) in spherical polar coordi- 
nates (r, 6, >) and using B = (BL, By, 0); Dungey [1954] 
obtained the following two equations relating the azimuthal 
component of the electric field of the disturbance EY to the 


azimuthal component of the fluid particle velocity V5 


ian 


(€) 

E | . 
sHi7 eee ae es EbOTaV RevTTA 10% rend iaivs 

i : 

tule ¥ (ares 
(Bh) (Fx Vv x ¥) e 


yitettmte. sd yam 1 bas tr war oe 


a haat ot 


anode ee.) seth Sas ai #4 aa 
lot Nvosds "6 aGTa¥ab oF beagnasie eee, i na 
bas 140 -$aer 240586, bre s ferigs aw 


ti 


a an “abe 
oh ut as Yi Pe 


nr a 


“fhto02 T6109 ea tsiiaz Hf (ey noi ssites eateemad . 
re cit ae ‘i ; ‘ 
[heen] yennud - (0 a) Pe) = s ghontey ts (e ee oss Na 


ere eta 
faddvmtse sat ontislar sho Fame on. wary il oni a Pets 
5 peste). es ! 
arid oF e sonnel ‘eds te bisit: he yo liad nen 


: e us foofay ar hine Shu 


42 


ene a° + 9° prea Ek 1M oa yr tsi. 0= E 
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Equation (5) is the equation for poloidal oscilla- 


tions (E,) which is coupled to toroidal oscillations (Vy) by 


*) 
the term on the right hand side of equation (5). Equation (6) 
is the equation for toroidal oscillations which is again 
coupled to poloidal oscillations by the term on the right 

hand side of equation (6). 

The above equations are too complicated to be of 
much use and they still have not been solved analytically. 
Dungey [1954] decoupled the two modes by assuming axial 
Symmetry (i.e., 5 = 0 which implies that the phenomenon is 
longitude independent and the disturbance may be considered 
to occur in phase over the whole earth). 


As a result of this decoupling procedure, the 


equation for poloidal oscillations becomes 
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(8) 
It can be shown from Ohm's Law and Faraday's Law 

and the axial symmetry consideration that the poloidal 
oscillations governed by equation (7) involve the following 


set of quantities 


and the toroidal oscillations involves the remaining set 


BEY Ey Gul eOR 0 b 4d; [0 0 Ved 


From the above sets of relationships, it can be 
seen that the plasma motion and the perturbation magnetic 
field lie in the meridian plane for the poloidal mode. The 
Poynting vector points across the field lines and the hydro- 
magnetic energy spreads to fill the volume. As a result, 
the whole cavity resonates and the period of oscillations 
for the poloidal mode will not be latitude-dependent. 

For the toroidal mode, the motion of plasma is in 


the azimuthal direction and the perturbation magnetic field 
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44 
is polarized in the east-west direction. The Poynting vector 
is directed along the field lines and the hydromagnetic 
energy is guided along them. Since equation (8) involves 
Operator (B,-v), the toroidal mode may be understood in 
terms of oscillations of magnetic lines of force where the 
period of oscillations will be latitude-dependent. 

In principle, the eigenperiods of the fundamental 
and higher modes of hydromagnetic oscillations may be computed > 
from equations (7) and (8). The first attempt to explain 
Pc's (under the old classification) in terms of poloidal 
oscillations was made by Kato and Akasofu [1955] who did not 
take into consideration the dipole character of the earth's 
magnetic field. This approach was further developed by Kato 
and Watanabe [1957] and Watanabe [1959]. These early analy- 
ses attempted to force the poloidal equation into a simply 
separable and soluble form through convenient assumptions 
concerning the plasma density and field dependence. The 
eigenperiods of oscillations depend upon the spherical shell 
bounded by the earth's surface and the outer boundary of the 
‘outer atmosphere’. A fundamental period of oscillation of 
about 180 seconds was estimated. The poloidal mode was also 
discussed later by other researchers such as Jacobs and 
Westphal [1964] and Carovillano et al. [1966]. The fundamen- 
tal period of poloidal oscillations is equal to four times 


the travel time for a hydromagnetic wave from the ionosphere 
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45 
to the boundary such as the magnetopause in the equatorial 


plane [Dungey and Southwood, 1970]. 

The fundamental period of the toroidal oscillation 
was computed by Dungey [1954] making use of equation (8). 
Using a constant plasma density distribution (p = 107! 8kg/em?) , 
and assuming perfect reflection at the earth's surface, he 


obtained the following relation for the fundamental periods 


where A is the geomagnetic latitude. T5 TOCA = 05. 1's) | 
minutes. Westphal and Jacobs [1962] derived similar equations 
for small amplitude hydromagnetic oscillations in cylindrical 
coordinates. Since the equations in this system are somewhat 
simpler than in the spherical coordinate system, they were 
able to compute the eigenperiods of toroidal oscillations for 
compressed dipole fields using a variable density distribu- 
mon. lf p varies as r~° where R is the radial distance from 
the center of the earth, vy, is constant and the eigenperiod 

of toroidal oscillation is equal to twice the travel time for 
a hydromagnetic wave moving along a field line between conju- 
gate points (as in the ‘stretched string’ model of Obayashi 
and Jacobs [1958].)..- For other forms:of plasma density distri- 
bution, the validity of this ‘time of flight® method in 


determining the fundamental period is reduced [Radoski, 1966]. 
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If the disturbance has a longitudinal dependence, 


it can be represented in the form ey where m is an integer. 
m = 0 corresponds to the symmetric modes described by the 
decoupled equations (7) and (8) presented earlier. Large m 
would correspond, of course, to the asymmetric modes which 
are associated with disturbances highly localized and varying 
rapidly in longitude. Dungey [1954] and Radoski [1967] have 
discussed the asymmetric poloidal mode. Since the perturba- 
tion magnetic field is transversely polarized radially in 
the meridian plane for this mode and its energy guided along 
field lines, it is referred to as the guided poloidal mode. 
The equations governing the guided poloidal mode and the 
toroidal mode are very similar. Radoski [1972] has considered 
the asymmetric toroidal mode and found that the periods of 
symmetric and asymmetric modes have a similar spatial depen- 
dence and are almost linearly related. Orr and Mathew [1971] 
have computed the eigenperiods for the symmetric toroidal and 
the asymmetric guided poloidal mode for dipole field lines. 
They found that the first harmonic period for the guided 
poloidal mode is approximately 20% longer than that for the 
toroidal mode and that higher harmonics have almost identical 
periods for the two modes. 

The effects of weak coupling between the symmetric 
modes through the inclusion of the Hall current in the 
generalized Ohm's Law and the use of a simplified field 


configuration (a hydromagnetic wedge model) have been con- 


ee mi 
| “aoabnaqed Ipaitbor tyno' Bat isinsdys er ong thn he 
Jsapedat ne atm ar ait arty wv} ant Cia reveal 
ont. uw bedi-yazeb iebenk a pang oii oF saab 
m sored TINE. boc nD2atg (8) Bae: 1%). snoiteupa i ‘git 
datiw 2ebow pivtaanges sft. ot 22799 to -onoqeoT7 06 
“priya has BestTedo! xlaein adnadtide th. aiFw bated: 
even [Vaer? thbober bees) [LP 284) vou mae abut tT pwol at rb 
+sdinired add sonte  . abot fébfofoa astydanamyes ont bee 
nt vlfarbes ben tra fog Ulaevavensad 27 a ere af = 
anofs babtue yorsns ett bis obon arnt Yor ene hg coth rem § 
otom febtatog bebiup sant <6 ot pahaten at i a 
sf bos Sion [shro Lag baht ue any oninnevon anorh 
boishteno> ena [SVOF)] rite obs v8 Limi 2) rev o8 nut: ta s 
| to ehoryeq sat ferris hque? bos abiom Pabbores ote 
-wogqeb leiigge vsTtire 6 avai 26b0m ohissamete, he 
[Iver] wansan bre 10 bets} 57 vive snh ly S2omte! 4x6, & : 92: 
bas fabrovd? 2 Pts ye arid 10% Paces vo, 
ead bist? fog b sot SBon 'ebtotag nabhug ‘in nemaye 
7 babiug sai “sto? boiva nossa dent oat dada dover 
sat 40% feds sn3 sa—noF 205 (let sintdougge zt abew Font >t 
Isottasbi Fzomis avert zat oumrnd saath Jang thee ~<a | 


a Y, Ss 2 alt 2abon ia ont vot 


hanes 


Di... 


atniemmye ens nssuiad ent Tewoo’ Asay Jo cists ont 
5 ahd ni dari (tah ate Boo Ae Ld dt se 
- tt on Paw 54 


Distt Beitiiqute 3 ©, oe 4 Fhe, nt) Bow 
Pee food ov si (1 sbom saben Saved ind a) no 


< ay 


47 
Sidered by McClay [1970]. He found that the toroidal mode 


resonance is spatially localized but the poloidal mode is 
also evident outside the resonance region. Because of 
coupling, a large portion of the cavity should oscillate 

with the same frequency. The poloidal mode exhibits close 
Spaced eigenfrequencies. The sense of polarization of the 
perturbations changes across the region of the toroidal 
resonance. Later, he studied the normal modes of oscilla- 
tions in a loaded concentric cylinder model [McClay, 1973]. 
In this model, the medium oscillates in its entirety without 
regard to local conditions and the eigenfrequencies depend 

on the arbitrary cavity size. Radoski [1974] has studied the 
asymptotic temporal behaviour of the asymmetric coupled modes 
in a cylindrical model. He found that as the system evolves 
toward the steady state, the energy is preferentially trans- 


ferred to the toroidal mode with the poloidal mode decaying 


away. 


2.3 Excitation Mechanisms 


In the previous section, we discussed the normal 
modes of oscillations in a hydromagnetic medium permeated by 
a non-uniform magnetic field when the system is perturbed. 
No mention was made of the sources of energy involved in the 
perturbation of the system. A discussion is therefore pre- 


sented here on the excitation mechanisms of long period 


micropulsations. 
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The most popular excitation mechanism is undoubted- 
ly the so-called Kelvin-Helmholtz instability, supposedly 
generated at the magnetopause due to interaction of the solar 
wind with the magnetosphere. The Kelvin-Helmholtz instability 
arises by the distortion of the interface between two fluids 
in relative motion. The first reference to the instability 
was made by Helmholtz [1868] with reference to hydrodynamic 
processes. Helmholtz's discussion was largely qualitative. 
The problem was treated analytically by Lord Kelvin in later 
years [Kelvin, 1910]. A simple hydrodynamic example of the 
instability is the phenomenon of ‘wind over water'. When the 
wind speed exceeds a critical limit (660 cm/sec according to 
Lord Kelvin), the instability will manifest itself as surface 
waves. Dungey [1954] was the first to suggest that the flow 
of the solar plasma along the boundary of the magnetosphere 
would generate waves in the same way that wind generates 
waves on water. The evidence from ground observations of 
micropulsations which suggest the Kelvin-Helmholtz instabi- 
ieyeas vaesource 15. the ditrwnalwanmiation pattern storm polar = 
zation in the horizontal plane (cc in the morning and cw in 
the evening). If surface waves are generated at the magneto- 
pause through the development of the Kelvin-Helmholtz insta- 
bility as the solar wind flows around the magnetosphere (to 
ve west on the morning side and to the east on the evening 


side as shown in Figure 12), then the plasma will have an 
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approximately elliptical motion, the rotation being in the 
Opposite sense along the dawn and dusk meridian. Since the 
magnetic field lines are 'frozen' into the plasma, the field 
lines near the surface of the magnetosphere will rotate 
with the plasma particles, thus generating elliptically 
polarized hydromagnetic waves which then propagate to the 
earth along the: magnetic field lines [Atkinson and Watanabe, 
1966]. The waves (looking along the direction of the field) 
will be polarized in the cc direction on the dawn side and 
cw On the evening side. This seems to explain the ground 
observations of diurnal variation of polarization in the 
horizontal plane, and the satellite observations made by 
Dungey and Southwood [1970], although there is some contro- 
versy on the stability of the magnetopause due to Kelvin- 
Helmholtz type [Sen, 1965; Dessler and Fejer, 1963; Fejer, 
1964; Southwood, 1968]. 

Swift [1967] has suggested that electrostatic 


oscillations may be responsible for long period micropulsa- 


tions. The electrostatic field couples to magnetic variations 


by inducing currents in the conducting ionosphere. The oscil- 
lation essentially occurs in a field-aligned column of 
enhanced (or depleted) ionization stretching between conjugate 


hemispheres. The energy comes from the free energy associated 


with gradients in particle density. 
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The interaction that occurs between hydromagnetic 
waves and energetic protons bouncing between hemispheres may 
lead to wave amplification [Dungey, 1965; Southwood et al., 
1969; Dungey and Southwood, 1970]. Quasitransverse modes of 
large azimuthal wave number m may be excited by taking energy 
from energetic particles through the bounce resonance 
instability. The change in energy W which results in wave 
growth is proportional to a change in the L shell parameter 
and is independent of the energy or pitch angle of the parti- 
cle. The energy exchange is a consequence of diffusion 
processes in (W, L) space. The most effective diffusion 
mechanism for the particles in L is the tilting of field 
lines in the meridian plane caused by the component of the 
wave field which is along the principal normal of the geo- 
magnetic field (e.g. the field of the guided poloidal wave). 
The mean velocity vector of the particle is also tilted. If 
the perturbation magnetic field reverses sign when the parti- 
Cheaounees, the particle zigzags across - as showngin Figure 
le. 

Kimura and Matsumoto [1968] suggested that Pc 5 
pulsations area result of the oscillations o1ethe siield 
lines due to a hydromagnetic instability which is excited by 
an electron or ion beam (such as precipitating auroral 
particles) passing through some part of the magnetospheric 


plasma. The excited wave has a frequency proportional to the 
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excess charge density, while the polarization of the excited 
wave depends largely on the sign of the excess charge. An 
electron rich beam causes instability in the Alfvén mode 
which is polarized in the cc direction looking in the irecmion 
of the magnetic field while a proton rich beam causes insta- 
Pity. in the modified Alfvén mode which is polarized in the. 
cw direction. However, there is a threshold beam velocity 
below which the instability can be suppressed. Kimura and 
Matsumoto claimed that the diurnal variation in the sense of 
polarization in the horizontal plane (i.e., cc in the morning 
and cw in the evening) can be interpreted in terms of these 
instabilities, although a mechanism in which different signs 
of excess charge will arise in the morning and evening was 
not discussed in detail by the authors. 

By extending the theory of the mirror instability 
[Chandrasekhar et al., 1958] to consider gradients in both 
the magnetic field B and the hot plasma density n, a co- 
existing cold plasma, and the effect of finite cyclotron 
radius, Hasegawa [1969] has suggested that the particle and 
field phenomena seen by Explorer 26 at L = 5 on April 18, 
E965 geomagnetic storm in the afternoon sector [Brown et al., 
1968] after the first large proton enhancements can be 
explained by the occurrence of a mirror instability in the 
magnetosphere. He has called this instability the drift 


mirror instability because of the coupling with the drift 
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Waves produced by VB and Vn.. For’ the onset of the insta- 
bility, it is necessary for the proton fluxes to have 
strongly anisotropic pitch-angle distributions and an 
enhanced energy density. When the perpendicular plasma 
pressure exceeds a critical value, the plasma tries to expel 
the magnetic field by its. diamagnetic effect. As a conse- 
quence, the plasma compensates for the decrease in the magne- 
tic field pressure and the parallel component of the local 
magnetic field is decreased. The instability has been called 
Pie mirnori dnsitability because:a lossi cone: distribution 
inherently creates such an anisotropic pressure. Since the 
mode causing the mirror instability is non-oscillatory, the 
oscillations are produced by coupling with the drift wave 
created by the ion drift perpendicular to the magnetic field. 
Thus, the drift mirror instability developed by Hasegawa 
explains the sudden inflection in the increase (decrease) of 
the proton fluxes (magnetic field intensity) and the subse- 
quent out-of-phase oscillations of the field and fluxes in 


the Pc 5 frequency range as observed by Explorer’ 26. 


Cladis [1971] has used a simplified model, based on 
a two-dimensional solution of Maxwell's equations, to demon- 
strate the possible occurrence of a resonance phenomenon in 
the magnetosphere. He has suggested that this resonance was 
responsible for certain correlated pulsations of the magnetic 


field and of the trapped particles observed in the outer magne- 
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tosphere (such as the one observed by Explorer 26 mentioned 
in the last paragraph). The resonance is excited by the VB 
drift of energetic trapped particles that have an inhomogen- 
eous distribution in the direction of the drift motion. The 
Field lines oscillate as standing waves in the meridional 
plane and form a wave pattern in the azimuthal direction that 
moves with a phase velocity equal to the mean azimuthal drift 
velocity of the particles. Particles with drift velocities 
near the velocity of this travelling wave become bunched in 
the potential wells of the wave. Particles with higher or 
lower drift velocities, and particles that drift in the 
Opposite direction, alternately lose and gain energy as they 
drift 'over' the wave. Consequently, these particles drift 
periodically across L shells, tending to follow the oscilla- 
ting field lines. Thus, the event observed by Explorer 26 
mentioned in the last paragraph can be explained within the 
context of Cladis' theory. The field line oscillations were 
initiated by the westward drift of the sharp proton front. 
The bunching of the low-energy protons accounts for the 
Osaiileations of the local. field. [he out-of-pnase, rejletion- 
ship of the proton flux and the magnetic field results from 


the diamagnetic effect of the bunched protons. 


2.4 Wave Coupling 


In section 2.2, we have discussed one group of 
theories of pulsations which deals with the resonance process 


and in section 2.3 we have discussed the other group of 
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theories of pulsations which concerns the excitation mechanism. 
A natural step to take is therefore to combine the active and 
passive aspects of the theories and treat the problem of 
coupling between the active and passive modes. Coupling 
between these modes has been studied in deter by Chen and 
Hasegawa [1974a,b], Hasegawa and Chen [1974], and Southwood 
[1974a,b]. 

Before we start discussing the coupling, some 
comments about the terminology of poloidal and toroidal 
modes are in order. Experimental evidence indicates that the 
assumption of axisymmetry which is used in many theoretical 
treatments is incorrect. For example, the assumption of 
axisymmetry is incapable of explaining the diurnal variation 
of polarization in the horizontal plane, the elliptical polar- 
ation and the spatially localized nature of the pulsations. 
Therefore, rather than assuming axisymmetry, the condition of 
| k 


are the wave numbers in the directions perpendicular and 


a| >> [ky | would be much more reasonable where ky and ky 
parallel to the magnetic field. For this condition, as 
pointed out by Lanzerotti and Fukunishi [1974], the distinc- 
tion between the terms poloidal and toroidal becomes meaning-— 
less because of the coupling of the compressional and shear 
Alfvén modes. It would therefore be better for discussions 


if the terms toroidal and poloidal were not used. 


The equation of motion of a fluid element is 


Ve ae ees 
peed Bo -evP (9) 
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where. ep is the mass density, V the perturbed fluid velocity, 


J the current density, B the magnetic flux density and P the 
plasma pressure. Assuming small perturbations, using 
Maxwell's equations V x B = Hod and substituting the dis- 
placement vector & (defined by S = V) for V, equation (9) 
becomes 


Chel saws b iene re 
Oe eer (vx b) x Bo + i (vx BL) x B - Vp (10) 
where b is the perturbation magnetic flux density, B the 
unperturbed magnetic flux density and p the perturbation 
plasma pressure. 
> > > 5 
Using b = V x (& x BY) and after some algebraic 


manipulations, equation (10) becomes 


Ve ra oz Bee i 
HP ane = (BY °V) 2 Ripe [p + US es Ga 


where 


€ = (Bev)B - (Bev) [(é*V)B + B(V-E)] 


Equation (11) indicates a coupling between a shear Alfvén 

wave (the right hand side is equal to zero) aud a surface wave 
whose dispersion relation is given by v(p + ood = 0. The 
strength of the coupling is decided by the vector € which 
appears as a consequence of the non-uniform magnetic field. 
Chen and Hasegawa [1974a] then solved the coupled equation in 


dipole coordinates while Hasegawa and Chen [1974] solved the 
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problem in Cartesian coordinates using straight magnetic 
field lines and obtained the essential feature of thee: 
results. 

In the simplified model of Hasegawa and Chen [1974], 
equation (11) becomes | 

ere yee Saeko east 
2 od UR SEO N ane ee MCSE GSE) 


a6 


(12) 


After relating the perturbations p and is through 
the adiabatic equation of state for a compressible fluid and 
the continuity equation, using b=V x (a x Be taking the 
Laplace time and Fourier space transformations of the dis- 
placement vector ig and assuming a field line resonance 
relatively localized in the east-west direction such that 
the azimuthal wavelength is much smaller than the length of 
the field line (i.e., ky >> ky), equation (12) further 


reduces to 


déé dine dé 2 
pe Yen = 
dy : dy aE by Y (13) 


Here the nonuniformity is taken in the y direction which 
corresponds to the radial direction directed toward the earth, 


the magnetic field is assumed to be parallel to the z- 
Cee 
direction and « = wu ely) - ky B (y) 


2 


Equation (13) describes surface waves (V by = #0) 


away from the resonant field line. Around the resonant field 
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59 
line, where € + 0, the second term in equation (13) dominates 
over the third term and a strong coupling between the surface 
wave and the resonant shear Alfvén wave occurs. 

Introducing a small positive ihddinars part of w 
ineensuth thatae = e, + ie., to take into account the 
ionospheric dissipation physically and to remove the singu- 


larity mathematically, [see also Southwood, 1974a] and 


expanding the equation around my ) = 0, equation (13) 


(@) 
becomes 
2 
doe dé 
—yi +—_1__7¥- ie = 0 (14) 
dy Soest Yea ey 
] (S35 
where y = ey (Southwood [1974a] has obtained a similar 
le 


equation for the perturbation electric field). 

The general solution of equation (14) can be 
expressed in terms of the moda etd Bessel function of the 
second kind. The corresponding eigenmode is shown to be 
continuous unless a monochromatic frequency source is 
assumed. An expression for E, can be obtained near the 
resonant point by using the incompressible and transverse 
nature of the shear Alfvén wave. The expressions for the 
horizontal components of the perturbation magnetic field B 
(which is related to & by 6 = ee BO)) can thus be 
derived and detailed predictions can be made concerning the 
sense of polarization and the orientation angle of the major 
axis of the polarization ellipse with respect to the point of 


resonance and azimuthal direction of wave propagation. 
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For negative azimuthal wave number (ky < 0) which 
corresponds to azimuthal propagation of the wave towards west 
from noon in the morning sector, the major axis lies in the 
second quadrant of the H-D plane if the Alfvén speed is 
decreasing toward the earth at the resonant field line. This 
tilt is predicted to switch across noon as the azimuthal 
propagation is towards the east from noon in the afternoon 
sector. 

The theories of Chen and Hasegawa [1974a] and 
Southwood [1974a] predict a reversal of sense of polarization 
in the H-D plane across the position of resonant field line 
and across noon as observed by Samson et al. [1971] (see 
Figure 8). 

Chen and Hasegawa [1974a] summarized EheeMn ve 
map tude, the ellipticity and sense of polarization as a 
function of radial position in the equatorial plane in a 
schematic diagram shown in Figure 14. A similar sketch was 
also derived by Southwood [1974a]. The diagram illustrates 
the essense of the theory of Chen and Hasegawa [1974a] and 
Southwood [1974a]: Surface waves (generated at the magneto- 
pause presumably by the Kelvin-Helmholtz instability) 
evanescing towards the earth are coupled to the shear Alfvén 
wave representing the oscillation of the resonant field line 


where the wave amplitude peaks sharply and across which the 


sense of polarization reverses. Hasegawa and Chen [1974] 
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have used the "guitar string" analogy to emphasize the 
monochromatic character of the surface wave. Picking one 
string of an infinite number of coupled Strings with a 
gradual change in lengths produces a wide band sound which 
corresponds to all the resonant frequencies of the strings. 
The sound decays as the inverse power of time due to phase 
mixing among oscillations of different strings. Hence, to 
excite a local field line, an excitor with a monochromatic 
frequency which corresponds to the resonant frequency of the 
particular field line is needed. 

However, a nonharmonic broad band source such as an 
impulse can excite sharply localized shear Alfvén waves at 
Sharp discontinuities in the plasma parameters (such as number 
density). This aspect of the theory for pulsations has been 
considered by Chen and Hasegawa [1974b] and Hasegawa and Chen 
[1974]. An initial. value approach was applied. They added a 
source term S which corresponds to the initial conditions on 


the right hand side of equation (13) which now becomes 


a 
abs d 
y Led lie Lekye Te 2 i 
dy” dye dy tf Cio ee (15) 


Chen and Hasegawa [1974b] noted that, as Uberoi 
[1972] had already pointed out, equation (15) has an exact 
analogy with the wave equation describing an electrostatic 
oscillation in a nonuniform cold plasma, which has been 


studied extensively by Barston [1964] and Sedlaéek [197la,b]. 
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63 
Using this analogy, Chen and Hasegawa [1974b] proposed the 
existence of a continuous frequency spectrum (non-collective 
mode) and a discrete frequency spectrum (collective mode). 

The non-collective mode corresponds to oscilla- 
tions with position-dependent frequencies and damping 
proportional to the inverse power of time. If the driving 
source has a continuous spectrum, oscillations with a 
frequency which varies continuously with latitude will be 
excited at all latitudes. Hence, to excite localized waves 
with a particular latitude-dependent frequency, a driving 
source composed of monochromatic waves is needed as men- 
tioned earlier for the steady state. 

The collective mode corresponds to the eigenmode 
of a surface wave excited at sharp gradients in the plasma 
distribution by an impulse which has a frequency spectrum 
that contains the eigenfrequency. The excited wave damps 
exponentially. The damping rate becomes smaller as the 
plasma discontinuity becomes sharper. As is expected from 
a surface wave, the collective mode has a peak in amplitude 
at the density gradient and decays away exponentially. As the 
profile becomes smoother and smoother, the mode is damped 
more and more until one is left eventually with a continuous 
spectrum. In the absence of a discontinuity, the eigenmode 
in the magnetosphere is simply a diffused damped mode 


corresponding to the non-collective oscillation and a mono- 


&3 ’ 
sity ‘Hazoqorg bannerd seve ei as ofens! 2 | 
‘sytasattoo- aon) muntosae vowsupayt) avon iano 6 ap 


. (sbon avi toetfoo) Ms DS: eneupen® Bist ay 


<BI Tigo Od ebpoge sT109 Sb dm ‘ey ivob hires oy edt 
Bi aoe ag pak spe TF ta inee nt 
oniwinb si FT amts +0 4swoe” aevount, eit Os Pei 
g nity 2noryet fraeo nega evounti nos ® eat 
5 ( if tw abuttdesl it iw clevounrtngs) @oPtay aati 
25Vew bashisdol -sitoks ad Spiek Jeobba test, ‘fs ree 
parvtab 6 von supsyt Snabasdeb ~sbud i tet sono ebony 
ast 3s hehbson 2f 2SV65W 5 tg emoytoomon 4a) bazognedig 
stehe vbsst2 eds Not “on 
sbomnopte end oF ebnogzertoo aivoin avisgal Téa. aate «! 
siesiq ont Ar 2inecbeye Grade 36, oe 
nurtasge CoOnsupeyt | rs 26 to Kat satuqnt ns yd 


eqmbb svew bottaxs ant arisiipe ri nee ts ene ented 

ont 26 val feme 2smaaad si6t i dy aT” wyirate 

mont bat aeqxe ey eke | Agaié te comoned witht inoaerB 
sbusrfams at 1590 5 zen bom ev bpefites ont .Svaiw. nie 
ait 2A A Est ansn0gxs, Vans, 218296 bis thathevw~ ie 


Bae Yentooua 2smo0ad. att 
ee 


beqmsb 2t sbom, ge i) che ie 


Zuounttnos F ad tw viteisqavs Fale 2% ‘ono fitev: siom & 
moet: PY ver erine es nae gis ey 

bubsepa Trib. b ania 2h aratqeoter 
-onom § bs tobielfeaso, “ev id09(tbd=nen 283 ot entont 


tT 


‘ak | 


oy 
me 


64 


chromatic oscillation is excited only by a monochromatic 
source. 

Chen and Hasegawa [1974b] have used the above 
theoretical approach to explain the observations of 
plasmapause-associated micropulsations observed by 
Lanzerotti et al. [1973] as shown in Figure 15. The damped 
oscillations observed at L = 3.2 are assumed to be the 
excitation of weakly damped surface eigenmodes which exist 
at the sharp density gradient of the plasmapause. At L = 4.0 
away from the plasmapause, only modes corresponding to the 
continuous spectrum were excited and, consequently, damped 
due to spatial phase mixing. Finally, it should be noted 
that damped Alfvén waves attributed to surface eigenmodes 
have also been observed in laboratory theta pinch experiments 


[Grossman and Tataronis, 1973]. 


2.5 Summary 


In this section, we shall evaluate the theories 
presented earlier in the light of relevant observational 
data. 

Since the time varying non-uniform magnetosphere is 
clearly not symmetric, the uncoupled symmetric poloidal and 
toroidal modes are not good approximations to actual condi- 
tions. As mentioned in the previous section, experimental 
evidence exists which demonstrates that the assumption of 


axisymmetry which is used in many theoretical treatments is 
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66 
incorrect. The symmetric modes are unable to explain the 
local time dependence of the observed pulsations, localiza-. 
tion in latitude and longitude of the pulsations, the 
diurnal variation in the sense of polarization, the ellip- 
tical polarization, the phase differences between stations, 
etc. As for the asymmetric modes, in order to maintain a 
mathematically tractable form models chosen for the study 
may not be particularly realistic. 

Although the ground observations of the sense of 
polarization in the horizontal plane in the morning sector 
and afternoon sector [e.g. Samson et al., 1971], the 
satellite observations of the local time dependence for the 
polarization near the magnetopause [Dungey and Southwood, 
1970] and the semiannual variation of geomagnetic activity 
[Boller and Stolov, 1970] seem to favour the Kelvin-Helmholtz 
instability, there is observational evidence against this 
mechanism as a source for micropulsations. Green [1976], using 
three stations at geomagnetic latitude  v 54° in the British 
Isles, found that the large majority of the Pc 3 and Pc 4 
pulsation events have the western stations leading in phase 
(independent of LT) which for a propagating disturbance 
implies eastward phase motion from the dawn side towards 
dusk. Such a direction of wave motion is contrary to 
Ghait! tof ta i erick wave generated by Kelvin-Helmholtz 
instability in the morning sector (i.e., westward 


propagation from noon towards dawn). Lanzerotti et al. 
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[1972] observed that three-quarters of the events detected at 


conjugate points near L = 4 during the winter solstice were 
cc polarized and that no significant switch in polarization 
was observed around local noon. These observations argued 
against Kelvin-Helmholtz instability as a source although it 
has been pointed out by Chen and Hasegawa [1974a] that the 
tilt of the major axis of the ellipse rather than the sense 
of rotation is a more critical parameter in finding the 
direction of azimuthal wave propagation at low latitude. 

Orr [1973] has also presented arguments against the Kelvin- 
Helmholtz instability as a cause for the dominant peak in 
occurrence frequency in the morning at a latitude of 70.9° 
observed by 01' [1963]. He has used the findings of Frank 
[1971], McDiarmid et al. [1972], and Alekseev and Shabansky 
[1972] regarding the latitudes of the last closed field lines 
to suggest that the pulsations reported by 0O1' were generated 
within the magnetosphere and not near the boundary so that 
the idea that the pulsations are due to surface waves at the 
boundary as a result of Kelvin-Helmholtz instability must be 
regarded as dubious. 

Swift [1967] explained the Pc 5 by an electrostatic 
drift wave instability. However, Hasegawa [1971] has pointed 
out that the electrostatic drift instability is rather 
unlikely to occur in the magnetosphere because of ion Landau 


damping in a high B plasma and the stabilizing effects of a 
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fractional mixture of cold electrons. In addition, the 
theory of Swift as it stands is developed purely for electro- 
static perturbations. Hence, a coupling scheme to the 
electromagnetic mode must be worked out to allow the theory 
to be evaluated using the observations. 

Green [1976] found that the longitudinal phase 
variation of Pc 3 and Pc 4 pulsations consists almost 
entirely of low m (azimuthal wave number) modes, which 
Suggest that the bounce resonance mechanism [Southwood et al., 
1969] is not a viable source. 

Although the theory based on the non-neutral plasma 
distribution [Kimura and Matsumoto, 1968] can explain the cc 
(in the horizontal plane) polarization in the morning and cw 
polarization in the afternoon, it is not apparent why differ- 
ent signs of excess charge should arise in the morning and 
evening. 

The local drift mirror instability [Hasegawa, 1969] 
may explain the storm time Pc 5 observed by satellites near 
the equatorial plane in the afternoon sector where the storm- 
time partial ring current is enhanced and the plasma pressure 
at the equator in the magnetosphere is sufficiently high to 
excite some of the high 6 plasma instabilities. On the other 


hand; Gladis [1971] has offered another plausible interpreta- 
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The steady state field line resonance model of 
Chen and Hasegawa [1974a] and Southwood [1974a] seems to 
explain the observations of Samson et al. [1971]. Figure 
8 which shows the diurnal and latitudinal variation of the 
Parizontal sense of polarization was compiled by Samson et 
oi Be [1971] using data from stations which were widely 
separated in the north. The two northern-most stations are 
separated by approximately 10° of latitude. This large 
separation between the two northern-most stations is crucial 
in determining the demarcation line and line of intensity 
maximum. Clearly, more stations are needed to fill in the 
gap in order to help distinguish among the many competing 
generation mechanisms for Pc 5 pulsations. Finally, it 
should be noted that Southwood [1975] has discussed the 
manner in which the field line resonance model may break 
down. He has pointed out the limitations of the linear theory 
which treats the continuously varying Alfvén speed across 
the field lines. 

The generation of surface waves at a density gra- 
dient by an impulse [Chen and Hasegawa, 1974b] explains the 
weakly damped oscillations observed at L = 3-2 by" Lamzerotti 
et al. [1973] and the damped type pulsations associated with 
sudden commencements and sudden impulses; i.e., Pc 4 and 
Psc 5 [Saito and Matsushita, 1967]. Thermal plasma density 
gradients exist at the plasmapause and in regions of detached 


plasma outside the plasmapause [Chappell], 1974]. The local 
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time distribution of regions of detached plasma is shown in 

Figure 16. It can be seen that the afternoon sector fing a 
‘large concentration of regions of detached plasma. Chappell 
[1974] pointed out that the smaller peak in the distribution 
of regions of detached plasma between 0800 andes EOOmL!. 1S not 
thought to be physically significant but probably represents 
an overlap in sampling in this local time sector. Therefore, 
the ‘cleaner! pulsations observed in the morning sector are 
probably not related to the surface eigenmode theory of Chen 
and Hasegawa [1974b] which requires a density gradient. 

Wilson [1966] has used a rotating line current as 
a model for the current systems driven in the ionosphere by 
hydromagnetic waves. However, it has been pointed out by 
Samson [1971] that this model does not explain the variation 
in the sense of polarization in the H-D plane with latitude 
and the large Z component at the latitude of the maximum in 
the intensity of the horizontal component. 

Finally, it should be noted that Hughes [1974] has 
studied the screening effect of the ionosphere and atmosphere 
on long period micropulsations. Assuming the ionosphere to 
be horizontally uniform, the flat earth to be perfectly 
conducting and neglecting the inhomogeneities in the magneto- 
sphere which cause coupling between the fast and the trans- 
verse mode, Hughes has obtained expressions relating the 


magnetic field of the hydromagnetic waves in the magneto- 
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Sphere to that observed on the ground. He has found that 
any disturbance seen on the ground due to a localized 
hydromagnetic wave in the magnetosphere will have no vertical 
current associated with it. Consequently, the magnetic field 
due to transverse waves must be Bach rotated through 90° 
when observed on earth, the rotation being caused by the 
efrects Of currents ‘flowing In the Tonosphere. If the 
horizontal scale length of variation is greater than the 
height of the ionosphere, the ‘polarization ellipse in the 
H-D plane observed at the ground must be rotated through 90° 
to obtain those in the magnetosphere. 

From the above considerations, it is not unreasona- 
ble to ade due that the theories for pulsations on the Pc 4 
and Pc 5 range are far from being complete. In this thesis, 
a new theory based on the temporal and spatial variations of 
the three dimensional current system flowing along the auroral 
oval is proposed. It will also be pointed out that the 
Steady state field line resonance model might be operative 


for some special and rare events. 
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CHAPTER III 
EXPERIMENTAL DETAIL AND DATA PROCESSING 


The data to be presented in this thesis were 
obtained from magnetometers operated by the University of 
Alberta along a common geomagnetic meridian. A line of 
these stations started operating in Western Canada along 
the corrected geomagnetic meridian ~ 301°E in the summer of 
1969. The line was in full operation in late 1971 and 
early 1972 with nine closely spaced stations ee rinantd the 
latitude range from ~ 60.6° (geomagnetic) to 83.1°N. The 
coordinates and code names of the nine stations are given in 
Table 2. The meieatian line of magnetometers discontinued 
Operations in the spring of 1972. In the Summer of 1974, 
four stations were set up in the region of the auroral oval; 
three of the stations were along a common geomagnetic lati- 
tude line. The coordinates and code names of the four 
Stations in operation during the summer of 1974 are given in 
Table 3. The locations of the stations are shown in Figure 
we 

Fach station was equipped with a three-component 
fluxgate magnetometer, together with a WWVB time code 
receiver, analog-to-digital system, and digital tape recorder. 


The detector head at all the stations except Resolute Bay 
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Code Names and Locations of the Magnetometer 


Stations in Operation During Late 1971 and Early 1972 


Station 


Resolute Bay 
Cambridge Bay 
Contwoyto Lake 
Fort Reliance 
FORt«Smith 
Fort Chipewyan 
Fort McMurray 
Meanook 

Leduc 


Station 


Code 


RESO 
CAMB 
CONT 
RELI 
SMIT 
FITCH 
MCMU 
MENK 
LEDU 
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was oriented in local magnetic (H,D,Z) coordinates. The 
detector head at Resolute Bay was oriented in the geographic 
(X,Y,Z) coordinates system. Data were recorded digitally on 
7-track magnetic tape at a rate of » 1 Sample/component/2 sec. 
The timing was accurate to within + 0.1 seconds. The dynamic 
range of the system was + 1000y with a sensitivity of they: 
The records were timed by recording WWVB directly on tape 
every 7 hr 38 min. A more complete description of the 
instrumentation is given in Appendix A. 

The data on the 7-track 'field tapes' were unpacked 
onto 9-track ‘master tapes' for each station. Where the data 
blocks could not be timed accurately by the computer due to 
the occurrence of bad data blocks between successive WWVB 
blocks, the. data were timed by correlating the micropulsations 
(in the D-component) occurring in this time interval with 
those recorded at adjacent stations. Much of the data recor- 
ded at Contwoyto Lake was timed manually in this manner with 
an caccuracy.of ~ + 1 minute. This particular process of 
timing proved to be a long and tedious task and phase infor- 
mation at the station of Contwoyto Lake 1s, of course,, less 
reliable. 

After the data from all stations were transferred 
to 9-track ‘master tapes', a series of ‘event tapes” were 
created with each tape containing a series of files. In each 


file all available data from the line of stations for a 
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specified period of time (event) were written using the same 
format as was used for the 'master tapes'. This allowed the 
Study of a specific event by only utilizing one file from an 
vent tape' rather than nine 'master tapes’. 

| The data in the event tape were then ready to be 
processed. The characteristics of the micropulsations con- 
Sidered in this thesis were determined by using the cross- 
Spectral matrix formalism. The basic concept of this approach 
was outlined by Born and Wolf [1959], and applied by Fowler 
ee al. £1967], Rankin and Kurtz [1970], Samson et al. [1971], 
and Fukunishi et al. [1975] among others. In our analysis, 
the data were detrended with a 1-20 mHz zero-phase shift 
Butterworth digital bandpass filter [Alpaslan, 1968]. A 
data window of length 60 minutes was utilized, and cosine 
tapers at either end of the data window were used to 
eliminate any remaining traces of the transient response of 
the filter and to improve the shape of the smoothed spectral 
window. The autopower and crosspower spectral estimates were 
calculated using the Fast Fourier Transform algorithm 
[Gentleman and Sande, 1966]. To employ the algorithm, zeroes 
were added so that the number of data points could be 
factored into integral products of 2. Raw and smoothed power 
spectra were plotted to determine the spectral peaks using 
the autopower spectral estimates. From the cross-spectral 
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79 
estimate at a frequency of component i and component j at 
the same station), the polarization ratios, the elliotica ties 
and polarization angles in the H-D, H-Z and D-Z planes were 
determined. From the station-to-station cross-spectral 
Matrices Oy ath cad (where a4 is the smoothed cross-spectral 
estimate for a frequency of component i at station if and 
component j at station aes coherencies and phases were 
determined. A more detailed discussion of the determination 
of polarization and phase from a three-dimensional vector 
time series is given in Appendix B. 

Since this thesis deals with the relationship 
between the micropulsation activity in the morning sector 
and the convection westward electrojet, information on the 
positions of the electrojet was necessary. To establish the 
latitudinal boundaries of the convection electrojet, the 
magnetic data were put in a latitude profile format. This 
presentation technique was used previously by Walker [1964], 
Bonnevier et al. [1970] and Kisabeth and Rostoker [1971]. A 
latitude profile shows the magnetic perturbations in three 
components plotted against latitude at a given instant of 
time. A quiet time baseline was chosen to determine the 
perturbation. The base line for each station was determined 
from the smoothed daily means of the magnetic variations 
which had excluded large magnetic perturbations. The lati- 


tude profiles are then interpreted using techniques developed 
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by Kisabeth [1972] to obtain information about the locations 
of the electrojet boundaries. A more detailed discussion of 
tne inferrence, of therdistribution ofselectric;currents in 
the ionosphere ane magnetosphere from ground magnetic data 
1s given in Appendix C. 

Finally, we note that the data were analysed in the 
local magnetic coordinate system (H, D, Z) for the polariza- 
tion characteristics of micropulsations and that the data 
were transformed to the centered dipole system (X, Y, Z) 
before being portrayed in the form of latitude profiles. 
Frequency bands which had distinct spectral peaks and 
high polarization ratios at most of the stations were 
selected for analysis. The intensity is the total power 
of the coherent signal (i.e., the sum of all three compo- 
ments°H; D andsZjectIif is the largest ae tWe intensities 
at all the stations in a given interval, then the relative 
intensity at a station i is defined by I,/1.. Intensity 
contours can then be plotted as a function of invariant 
latitude and Universal Time. Universal Time is used through- 
out the text. An approximate conversion to local geomagnetic 
time at our stations can be made by subtracting 85 hours from 
the given Universal Time while local time is approximately 
Universal Time less 7 hours. The date of an event is 


expressed in terms of number of days after the beginning of 


the year (January 1 is Day 1). 
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CHAPTER IV 
Pe 5 ACTIVITY IN THE MORNING SECTOR 


In this chapter, we fired present observations that 
the micropulsation activity in the Pc 5 range in the morning 
sector is closely related to the convection westward elec- 
trojet. The evidence that the electrojet itself exhibits 
Spatial oscillations are dealt with next. From the observa- 
tional characteristics of Pc 5 micropulsation activity as 
related to the westward electrojet, a model regarding the 


generation mechanism of the Pc 5 activity is then presented. 


4.1 The Relationship of Pc 5 Activity in the Morning 
Sector to the Auroral Westward Electrojet 

This section includes and extends the work of 
Lam and Rostoker [1975] reported earlier. The results of 
analysing three days' activities will be presented to 
meuStrate the relationship of the Pc Ss activity “in. the 
morning sector to the auroral westward electrojet. 

The events selected for analysis are chosen 
from three separate days. They are Day 253, 1971 and Days 16 
aneatl?y.-1972.. All the events occurred inthe morning sector, 
spanning the local time regime from pre-dawn to local 
The magnetograms both unfiltered and filtered 


magnetic noon. 


(with a 1-20 mHz digital filter) for the events are shown in 
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82 
micures: 185. 19, 20 (Day-235, 1971), 21, 22 (Day-16, 1972). 
23 and 24 (Day 17, 1972). A magnetogram is an amplitude- 
time plot of the magnetic perturbations in the H (the 
component pointing towards the local magnetic north), D 
(the component pointing eastward) and Z (the component 
pointing downward) coordinates recorded at our Stations, 
the northern-most ‘station being RESO (see Table 2 for the 
name of the station) and the most southern one being LEDU. 
The perturbations are plotted as a function of universal 
time (local time for our stations is approximately equal to 
universal time less seven hours). It is apparent from the 
magnetograms that the amplitudes of the oscillations are 
rather large (up to ~ 100y peak to peak for Days 16 andi17 ). 
The pulsations seem to occur in bursts and over hit tenent 
latitudinal ranges at different intervals of time. For 
example, in Figure 24, one can see a northward shift of 
pctayity asia stunction oftimeg Jhis.trendiwill later be 
Shown to be related to the locations of the convection 
westward electrojet. 

The locations of the convection westward electrojet 
plotted as a function of invariant latitudes and universal 
time’ will be used as a frame’ of reference for\tne presentation 
Of tthe) pulsation data. Both the poleward border and the 
equatorward border of the electrojet were inferred from a 
series of successive latitude profiles of the magnetic 


perturbations (see Appendix C for more detail regarding the 
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mnterpretations of the latitude profiles). 

Some examples of the latitude profiles from 
which the borders of the electrojet were inferred are shown 
iietigure.25. . the, borders of the electrojet are indicated 
by arrowheads in Figure 25 and were inferred by assuming 
that there is a gaussian distribution current density across 
the electrojet. For each profile the horizontal components 
(H,D) were tee to provide orthogonal horizontal compo- 
nents (X,Y) in the centered dipole system. For our stations 
dipole latitude is essentially the same as invariant latitude 
(invariant latitude was chosen because the project was 
initiated with a view of comparing ground pulsations data 
with satellite data). 

Figure 25A shows that the convection westward 
electrojet was located near our southern stations around 
Poy UT. About san hourslater at 1631 UT this.wWestward 
electrojet was located over our northern stations as shown 
Poeearoures25B% Figure.25C and Diillustrate {the vapid 
poleward motion of the equatorward border of the westward 
electrojet within ayshort time. slit can be seen that the 
equatorward border of the westward electrojet shifted 
northward substantially (about 4°) in a few minutes. It is 
interesting to note that an eastward electrojet (Figure 25D) 
appeared in the latitudes where the lower portion of the 
westward electrojet was (Figure 25C). Finally it should be 


noted from these latitude profiles that the Y profiles which 
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exhibit a more or less step-like function (Figure 25B shows 
the step clearly) are indicative of the net downward field- 
aligned current flowing into the latitudinal regime 
occupied by the westward electrojet. The Cie of ine! 
electrojet with respect to our line of stations probably 
prevents profiles in Figure 25A,C and D from showing a clear 
step-like function jin Y. 

To see how.the Pc 5 activity is related to. the 
location of the westward electrojet, intensity contours of 
a selected frequency band were plotted against the back- 
ground of the electrojet as a function of invariant latitude 
and universal time as shown in the next few figures. Each 
figure consists of two portions. The upper portion is the 
intensity contours of the micropulsation activity plotted as 
a function of latitude and universal time. The stars on 
tne ordinate indicate the locations of our stations in 
invariant batitude. Only the first andMleast tetrver of the 
Goce: fname ofsour stations<are yplotted. Jror example, RO 
represents RESO (Resolute Bay)» The lower*portion of. the 
figure consists of the poleward borders and equatorward 
borders represented by solid lines of the electrojet (which 
were inferred from a series of latitude profiles about 10 
minutes apart) and the intensity maximum of the pulsation 
activity represented by the shaded area. 

Figure 26 shows the intensity contours of the 


pulsational activity in the lower frequency band CTs thae 5 
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mzjeoft.the Pc 5 spectrum -for Day 16,1972. The 'meandering' 
of both the electrojet and the intensity maximum with time 
can clearly be seen. The next two diagrams, Figures 27 

and 28, illustrate that the Pc 5 activity in the higher 
meequency bands of 2.0 - 3.0 mHz and 3.3.- 5.1 mHz. ‘for 

Day 235, 1971 and Day 17, 1971 respectively again " follows’ 
the electrojet very well. The proximities of the contour 
lines in the upper portion of the above figures and the 
locations of the intensity maximum with respect to the 

| electrojet Shown in the lower portion of the figures for 
different Spectra bands in the Pc 5 range for three 
different events clearly illustrate that the region of Pc 5 
oscillations is confined to the region of the auroral 
eliectrojet. 

To demonstrate further that the micropulsation 
activity is intimately related to the westward electrojet, 
we present in the next few figures the power spectra 
corresponding to different locations of the westward 
electrojet. The powers in H, D.and Z covering a frequency 
range of 1 to 16 mHz were plotted for all our stations, 
RESO being our most northern station and LEDU our most 
southern station. The power spectra are characterized by 
six degrees of freedom. The spectrum is truncated at 
Oa yo /mHz power level because of possible quantizing errors 
Figure 29 shows the power spectra 


at lower power levels. 


for the interval 1100 - 1200 UT for Day 17, 1972. During 
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this interval, the westward electrojet was located near 

our southern stations. The stations that were within the 
electrojet most of the time are indicated in the figure. 
Mevcan be iseen that inside the electrojet and near the 
equatorward border the power level was much enhanced, 
particularly in the H component. The presence of high 
frequency peaks in the Pc 4 band in the stations inside the 
electrojet was evident. The stations poleward of the 
electrojet have very low power and show no high frequency 
peaks at all. However, as the electrojet shifted poleward, 
the enhanced power levels and the high frequency peaks also 
Shifted northward as shown in Figures 30, 31 and 32. The 
contrast can clearly be seen in the power spectra for CONT 
for the four different time intervals. 

To demonstrate further that the pulsational 
activity peaks in the same latitudinal range occupied by the 
westward electrojet and to show that all the spectral 
components in the Pc 5 band peak in the same latitudinal 
range, we present in the next three figures the latitude 
protiles*ery the powers of different frequencies. for three 
time intervals corresponding to three different locations of 
the westward electrojet. The magnitudes of the powers are 
not shown in the diagrams because we are only interested in 


where the powers of a particular frequency band peak. Be- 


sides, the powers almost drop to noise, levels outside the 
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confines of the electrojet as is evident in the intensity 
contour plots (see for example Figure 28). Nonetheless, 
the magnitudes of the peak powers were compiled in Table 4 
which gives an idea of the relative intensity of activity 
in the different frequencies and a comparison of the 
magnitude of powers in different components. 

It can be seen in Figure 33 that when the westward 
electrojet was located at low latitudes (the dashed lines 
indicate the borders of the westward electrojet), all the 
Spectral components (1.7 mHz, 5.0 mHz and 7.1 mHz) peak 
within the confines of the electrojet. As the electrojet 
Shifted to higher latitudes as shown in Figures 34 and 35, 
the peaks of powers in different frequency bands also 
Shifted. Again note the contrast in power level at CONT 
(the third dot from the right) for these three time inter- 
vals. These observations that the low frequency end of the 
Pc 5 spectrum maximizes at low latitudes (~ 64°) as well as 
apenigh. batitudes (~ 73°) and that differentuspectrat 
components all peak in the same latitudinal range occupied 
by the westward electrojet clearly suggests that the current 
System associated with the convection auroral electrojet 
dominates other sources (such as the resonant oscillations 


of field lines) in its contribution to the low frequency 


band of the micropulsation spectrum. 
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The close relationship between the Pc 5 activity 
and the convection westward electrojet can further be 
substantiated by the results of the polarization analysis 
which now follows. The shapes and tilts of the polarization 
Peli ipses as well as the sense of polarization (see Appendix 
B for more details regarding the determination of the 
polarization parameters from a three dimensional vector time 
series) were projected on to the three orthogonal planes. 
The plane is viewed from a direction in which the 
coordinates form a left handed system (e.g., downward on 
oe horizontal H-D plane). The polarization ellipses in the 
H-D, H-Z and D-Z planes for the three different spectral 
bands chosen from the three events we presented were 
positioned as a function of latitude and universal time with 
the convection westward electrojet Superimposed on the back- 
Ground as in Figures 36, 37, 38 (Day 16, 1972 in the 1.3-2.5 
waepand), 39, 40, 41 (Day 235,_19714in the 2.0-3:0-mHz band), 
WeeeAs and 44 (Day 17, 1972 in the 3.3-5.1 mHz band). ~The 
eumapses ‘that are shaded,indicate polarization Jn the ilock- 
wise direction (CW) in the plane while those that are not 
shaded are polarized in the counter-clockwise (CC) direction. 

An examination of the polarization ellipses in the 
H-D plane (Figures 36, 39, and 42) would indicate that tne 
behaviour of the ellipses follows the electrojet. There 


appears to be a tendency for the ellipses to be tilted in the 
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second quadrant of the H-D plane inside and in the vicinity 

of the westward electrojet, indicating that the phase differ- 
ence™between H and D®at’a particular locationcinorrclose to 
the electrojet is greater than 90°. Away from the electrojet 
in the south, the ellipses have a tendency to be tilted in the 
first quadrant of the H-D plane indicating that tne phase 
difference between H and D is less than 90°. The ellipses 

are more inclined towards H than towards D implying that 

the perturbations in H are generally greater than those in D 
for the Pc 5 micropulsations. The sense of polarization for 
the=stations in the south is genérally CC. (indicatingethat 

H lags D in phase) while the sense of polarization for the 
stations in the north is predominantly CW (indicating that H 
leads D in’ phase). This pattern agrees with the polarization 
analyses by previous workers (see Figure 8). However, with 
the locations of the convection westward electrojet as a 

frame of reference, our analyses indicate that the demarcation 
line (the line across which the sense of polarization rever- 
ses) follows the electrojet, shifting northward with the 
electrojet (Figures 39 and 42) or ‘meandering’ with the 
electrojet (Figure 36). Note that in Figure 36, for the UT 
Phitey vale 1500-1600,"the polarization’eblipses at the southern 
Stations have CW sense of polarization while’ thersense of 


polarization at the northern stathonsmirs? CClLeerhise particular 


pattern appears to be opposite to the statistically determined 
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Pomarizdtion pattern shown in Figure 8hfor:thisopanticular 
local time. However, this happens at a time, when the 
westward electrojet has reached its furthest southward 
position. Based on previous results, one would anticipate 
the senses of polarization. to switch to CC equatorward of 
BevULSoSthat the switch in senserof polarization ine the. H-D 
plane would occur at the equatorward border of the electrojet. 

In the H-Z plane (Figures 37, 40, 43), the sense 
Of Dolarization of the ellipsesnseemstro be opposite to 
that in the H-D plane with predominantly CW polarization in 
the south and CC polarization in the north. Again, the 
demarcation line appears to follow the electrojet. The 
Eaiapses inside and in the vicinity of sthe electrogetjare 
mostly tilted in the second quadrant indicating that the 
phase difference between H and Z at a station is greater 
Phan 90° whilethe: eltvipses far awaye from thepelectroget, are 
tilted predominantly in the first Guadrant indi catinige tivat 
the phase difference between H and Z is less than 90°. The 
ellipses are inclined more towards H than towards Z implying 
that the H perturbations are greater than the Z perturbations. 
H lags Z in phase for the CC poliar izeti one ell sese andsclieads 
Z in phase for the CW polarization eld ipsiess 

The polarization in the D-Z plane (Figures! 38,41, 
and 44) appears to be a little more comp ipeaivedi.7a Wewe nthe less. 
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tmere appear to be two reversals in the sense of polariza- 
tion, one associated with the poleward border of tne electro- 
jet and the other associated with the equatorward border. 
Going along our line of stations, one encounters CW ellipses 
(D leads £if-pnase)y. then; CO ellipses (D tags.2 in) phase) 
anid CW ellipses again. In general, the ellipses seem to be 
inclined at an angle of 45° towards the axes indicating that 
the magnitudes of perturbation in D are approximately equal 
meenose in Zz. It seems that the tilts of ellipses inside 
the electrojet are sone to those outside of the electro- 
foeoeri gure 44 illustrates this statement best). ‘The tilt 
in the first quadrant implies that tne phase difference 
between D and Z at a particular station is less than 90° 
while the tilt in the second quadrant implies that the pnase 
difference between the two components is greater than 90°. 

We summarize tne general senses of polarization 
Senematically using tne results of Day 1/.,\1972 in» the 
following three diagrams. 

Einoure. 45 shows the “sense of polarization in the 
feoepbane. -The point-that this diagram is intended to 
iulustrate is. that the demarcation line follows: the electro- 
jet. 

Figure 46 shows the sense of polarization in the 
H-Z plane and indicates that the sense of polarization in 


the H-Z plane seems to be opposite to that in the H-D plane. 
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Figure 47 shows the sense of polarization in the 
D-Z plane. Except for the time intervals when the wnole 
eueetrojet undergoes’ rapid variations in-its position (for 
example the time interval of 1500-1700 on Day 16, 1972 shown 
in Figure 38), there are two reversals in the sense of 
polarization, one associated witn the poleward border of the 
electrojet, the other with the equatorward border. 

The senses of polarization discussed above can be 
understood in terms of the relative phase changes across the 
Stations. Typical examples of the phase changes are shown 
in the next two diagrams. The top three panels in Figures 
48 and 49 show that the powers in three components peak within 
the confines of the electrojet as mentioned earlier. The 
bottom three sere show the relative phase changes across 
tie Sitations in three components. “The phase difference 
between any two components at a given station is not snown, 
aseainrorhnation on this, 1S.impltiert in thessiapes se Cin tsand 
Sense of polarization of the polarization ellipses Shown 
earlier. The magnitude of the phase at a given station is not 
important because of its cyclic nature. teas eeneronot ile OF 
the phase that is of interest. It cones sSeenai rom Figures 
48 and 49 that there is one rapid phase change of about 180° 
in H across the region of maximum power; there is very little 
phase change in D; there are, two rapid phase changes in Z. 


Thus, the reversal in the sense of polarization in the H-D 
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plane can be attributed to the change in phase of tl and the 

mwooreversals in the D-% plane can readily be understood as 

eter to phase changes’ in Z.. The sense of polarizatton in. the 
H-Z plane can be understood using information: regarding -the 

phase difference between H and Z at the stations in conjunc- 
tion with the eleriwe phase cnanges across the stations. 

We have shown earlier in tnis section feen the 
intensity contour plots that the intensity maximum follows 
the electrojet. The results of the polarization analysis 
eeemindicate that the demarcation line follows the electro- 
jet. It would therefore be of interest to see now the line 
of peak intensity is related to the demarcation line particu- 
Tariyein ithe horizontal plane. The’ results of superimposing 
theshine of peak intensity and the demarcation» line in H-D 
pane; on the background electrojetyfon the tnree events=are 
Soom in Figures. 50, 51 and/52. Its apparent trom these 
pilots, that, while both lines follow,each, other; they do nor 
coincide with each otner, which implies that linear polariza- 
tion does not occur in the region of peak intensity. We feel 
that our stations were sufficiently close enough to render the 
above statement valid. 

We conclude this section by showing how pulsations 
respond to electrojet motion. We have compiled in Table 5 


the hourly averages of interplanetary magnetic field parame- 


ters and the Kp values. It can be seen from Table 5 and the 
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MABELE S~ Hourly averages of interplanetary magnetic field 
parameters (6 is the angle between the field vector 
andthe solar ecliptic plane im dearee, and B. is the 
north-south component of the field in y) and Kp values 
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plots containing the positions of the electrojet that the 
electrojet moved in response to the changes in the inter- 
planetary magnetic field reflecting the changing level of 
magnetospheric convection. The Kp values indicate that 
‘the Magnetic conditions for our events were moderate. 

: Figure 53 shows that the pulsational activity in 
mene ¢c-3 mz band jumped from low signal levels to peak 
intensity when the electrojet was in a transient state 
moving from low latitudes to higher latitudes. A second 
maximum, though not as intense as the first one, occurs after 
the electrojet reached its northern positions. The inter- 
planetary data indicates that there were changes in the 
interplanetary magnetic field in the time interval for the 
second maximum (viz., the normal component of the IMF turned 
southward in the interval 16-17 UT). Thus, although the spa- 
tial condition of the electrojet was stable, the change in 
level of magnetospheric convection due to changes in the 
interplanetary magnetic field might produce changes in the 
imgensity of the current system-associated withy thesetectro- 


get resulting in the burst of activity. 


Figures 54 and 55 again show the enhancement of 
signal strength in the low frequency band and higher 
frequency band of the Pc 5 spectrum associated with changes 
in electrojet configuration. The higher frequency band 


appears to respond earlier to motions of the electrojet than 


the lower frequency band. 
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Figures’ 56,and 57 11lustrate that both. thes how 


frequency band and high frequency band of the activity on 
Day 1/7, 1.972 jumped to 20%, of -their,peak intensities as the 
equatorward border moved northward rapidly. The signal 
evel then remained relatively Unchanged over the next three 
hours until the time when there was another rapid motion of 
the equatorward border northward. The pulsational activity 


reaches its peak power at that time. 


A°2 Evidence for the Spatial Oscillations of the Electrojet 
In the previous section we have presented the 
characteristics of pulsations with the electrojet super- 
imposed on the background. The locations of the electrojet 
boundaries were inferred from a series of successive latitude 
profiles at about 10 minute intervals. The result of plotting 
the locations of the borders as a function of time would 
therefore only indicate the trend of the electrojet giving 
only a macroscopic picture of the Wocdtions to7 > tiemeniectro— 
jet. We now present in this section the locations of the 
electrojet which were inferred from a series of successive 
latitude profiles only 23 seconds, apart. Thus the result 
of plotting the borders of the electrojet as a function of 
time would provide a microscopic picture of the spatial 


behaviour of the electrojet. It will be shown in the 


following diagrams that the borders of the electrojet behaved 


an an oscillatory fashion and that the os ci biazions, cof tne 
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parders seem to occur in unison with the micropulsations.. 

The magnetogram shown in Figure 23 for Day 17, 
1972 event shows that in the interval of 1300 and 1400 UT, 
the pulsations a high latitude notably at CONT are of Taare 
frequency than the pulsations at lower latitudes. This is 
due to the fact that the poleward border of the electrojet 
oscillates more slowly than the equatorward border of the 
prectrojet as illustrated in-Figure 58. In Figure 58, we 
have plotted the poleward border and the equatorward border 
(inferred fnom latitude profiles 23 seconds apart) as a 
function of Universal Time. We have also superposed the 
pulsations in the Z component recorded at CONT (which is 
close to the poleward border) and at FITCH (which is near 
the equatorward border). It can be seen that there is a 
remarkable resemblance between the border oscillations and 
magnetic oscillations. As mentioned earlier, the timing for 
CONT might be off by as much as 1 minute. However, this error 
on the very long period wave trains shown in the diagram would 


produce effects which are hardly visible to the eye. 


Figure 59 shows the correspondence between the 
motion of the equatorward border and the pulsations recorded 
at MCMU and MENK just slightly south of that border. Note 
the higher frequency oscillations recorded at v 1220 UT seem 
to correspond to the higher frequency oscillation of the 


equatorward border of the electrojet. 


Figure 60 illustrates that both the poleward and 
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equatorward borders oscillate in a similar manner for the 
interval 1700-1800 UT in contrast to the interval 1300-1400 
UT. As a result, the station close to the poleward border 
(ies, CONT) and the station close to the equatorward border 
(i.e., RELI) oscillate with the same frequency in contrast 
to the oscillations recorded at CONT and FTCH between 1300 ~ 
and 1400 UT which are quite different. 

The border of the electrojet which is determined 
from the AZ extremum from the latitude profile seems to 
oscillate in phase with the magnetic pulsations in the Z 
component and slightly out of phase with the horizontal 
components particularly with the D component. This is 
illustrated in Figure 61 where the magnetic pulsations in 
3 components recorded at SMIT are presented together with 
the variations in the position of the equatorward border of 
the electrojet. 

Avplot’ oT. the magnitude of tonetnegative Ax extre- 
mum (denoted by HMAX) of the latitude profile as a function 
of Universal Time is shown at the bottom of Figure 61. It 
indicates that the three dimensional current system producing 
the H perturbation also fluctuates in: Tts. intensity. | hus’ it 
would appear that the Pc 5 micropulsations contain contribu- 
tions from both the spatial oscillations of the westward 
electrojet and temporal oscillations in the intensity of the 


electrojet. In the following section we Shall explore the 


implications of this suggestion. 
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4.3 Generation Mechanism of Pc 5 Micropulsations 

Our study of the latitudinal variations of the 
spectral components of the Pc 5 micropulsations indicates 
that all spectral components peak at the same latitudinal 
range and that the low frequency component maximizes both 
at low latitude and high latitude (see Figures 33, 34 and 
35). It is evident that these observations do not readily 
Support the field line resonance model which predicts a 
latitudinal dependence of the spectral components. Further- 
more, it has also been shown that the demarcation line in 
the horizontal plane does not coincide with the line of 
intensity maximum as predicted by the field line resonance 
model (see Figures 14, 50, 51 and 52). Besides, there has 
been no conclusive experimental proof for the existence of 
Kelvin Helmholtz instability which has long been regarded 
dasa source ifor .the .excitation of fiel dhki nepnesonanceyamlle 
are therefore compelled to consider other alternatives for 
the generation of Pc 5 micropulsations although it must be 
emphasized that the field line resonance model might be 
operative in some specific occasions as will be shown in the 
next .chapter. 

In the previous sections, we have presented the 


pulsation data against the background of the convection 


westward electrojet. It is obvious that the intensity 


maximum and the polarization follow the electrojet (see for 
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example, Figures 26 and 47) and that the Signal strength 
of the pulsations is markedly enhanced in conjunction with 
sudden changes in the position of the equatorward border of 
the electrojet (see for example, Figure 57). We have further 
demonstrated that during and after the sudden shifts in 
position, the equatorward border of the electrojet oscillates 
spatially at a frequency which is a dominant component of 
the Pc 5 spectrum (see for example, Figure 59). These facts, 
coupled with the observation that there is always a large 
perturbation,in Z, point to the convection westward electro- 
jet and its associated three dimensional current system as a 
possible source for Pc 5 pulsations. 

fo obtain a rough picture of the physics¥of the 
current system we may regard it as an electric cireuit with 
lumped impedances as shown in Figures62.-"Bostroma.l 972) aes 
used a similar equivalent circuit for the substorm current 
system and found that the circuit rings sin the shinee frequency 
range in response to variations in current flow through the 
jlonosphere. The ringing mentioned above represents LC- 
oscillations of the equivalent circuit, and we shall follow 
Bostrém's technique by attempting to construct an equivalent 
circuit whose ringing will constitute Pc 5 micropulsations. 
Holzer and Reid [1975] and Reid and Holzer [1975] have also 
used a similar circuit analogy to study the response of the 
coupled magnetosphere-ionosphere system to changes in dayside 
field line reconnection rate. The construction of our equiva- 


lent circuit is based on the following reasoning. 


ddpnsite [pnp at 


ddfw no cae a nt ‘boo nnsins wit boy 


mt > bis os | 

estehi roe Stoxd sets: eid a0 na | 
+6 Tnonoqnes sngatmob Py et ‘Hote. yomaupert 6 te N 
,2s2AF S29AT/ ol Re up Ta. {quexs et 982) ree , 
op yet 5 ll i Sang, teas apap te 
| Sintog <S' nf Mobiet 

5 25 MOS 2X2 inserts) | cnot enamtb Se4d bersiz0ce0 ww 
| 290} tae tag e a 0% sowie 6 rd 
afd Vo 2sreyng ait to auaate ‘tower B aketdo! opwed 

(iW dtusvts optgoe Es ns 20 sae Vea ‘ow wing eR bhi 


A ; m nb e 
oe my foe, iM aa a ne et a 
26h £SVOTP mon $ 2 38 Bei: 
i inte A : 


i ei th nwo eB seen 


-ortasts by swh2ew eotiaavne9 ait o 


ies 


Yonsiipsyt § Henly 


eng ipuorda \ 


3 ednge ae d : oo nom 
gee Ee i! ele : 
wolfot Thad, ow. ne a iil piven ant Fo | 


14, 


tusfevi tups 6 “fountantaa os | cent 


eno tt s2/uqoxa Ft a af at 09 Hw atent perl pi 4 
oats ayved raver] vaso bas sf ns ; fever]: bref bas! 49sTot 4 


Si. oe teotTs) sinha 
} rn z 2 any oityzot 


e.. uP re ie 


-sVFUpS tho to not d2ix3 ena Ligh i am Bet? senn cae oni oe 


sj to eznogesy oad bude oF 


aht2yeb ni eepnada of datevee 


R 


Figure 62 


» oo 4 
5 an 
AM ke ele) ie 
; mera 
n : 
ihe ; 
‘, 
f 
: 
; 
® 
i ~4 
am 
‘ 


oy : mu va iy Aci 
ab) ete 
4) agi Hy by ‘nap 
er ee 
i ‘' ; , Moe 0 4 L. 
hy Las W 1 ae 
mA Lan: x 


147 


In Figure 5 we have shown the configuration of 
the auroral electrojets and field aligned currents associa- 
ted with convective processes in the magnetosphere. From 
this figure it would appear that the westward electrojet 
stretches from near noon to midnight. The huge scale size 
of this system makes it implausible that Pc 5 micropulsations 
could be explained by LC-oscillations of any electrical 
circuit associated with the volume of space occupied by the 
three dimensional current system. For example, Rostoker and 
Bostrom [1975] calculated values of L = 35 H and C = 9 x 107° F 
for a smaller circuit associated with the eastward electro- 
jet in the evening sector and found the frequency of LC- 
oscillations to be ~ 0.1 mHz which would be well outside 
the Pc 5 frequency range. However, Hughes and Rostoker 
[1976] have pointed out that the ionospheric conductivity 
discontinuity at the dawn meridian must be associated with 
the divergerce of desiward current in the ionosphere into 
field aligned current flowing up the tield lines iN corn une 
magnetosphere. Thus there is, confined to the pre-noon 
Guaardant., a three dimensional current system involving 
downward field-aligned current near noon diverging into 
westward ionospheric Hall current with subsequent CV en 
into upward field aligned current flow at the dawn 


gence 


meridian. It is oscillations of this current system which 


we shall use to attempt to explain Che Orrvgr, ot Fre eo 


micropulsations. 
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Consider the dawn and dusk meridian and the con- _ 
vecting plasma drifting through it. Rostoker and Bostrom 
[1976] have proposed a mechanism by which the Birkeland 
currents are driven through the extraction of kinetic energy 
from the magnetotail plasma by braking the motion of the 
plasma as it drifts towards the flanks of the magnetosphere. 
Moder and Lucht [1974] have found that the equatorward 
directed electric field decreases in magnitude as One moves 
from pre-dawn to post-noon (see Figure 63). Hence, as the 
plasma moves from a tail like to a dipole like magnetic field 
configuration as shown in Figure 64, it is slowed down as 
indicated by the decrease in the equatorward directed 
electric field. One may then consider that the energy 
extracted from the drifting plasma may provide the electri- 
cal Be 6 drive current systems in the pre-noon quadrant. 
Hence, the voltage source in Figure 62 is associated with 
the extraction of energy from the convecting plasma. Con- 
sider a region of space bounded by field lines poleward and 
equatorward of the auroral oval and confined within three 
time zones. Then, the current carrying volume can be 
represented by a rectangular box of length & (length of field 
lines in the generator region) ~ 15 Rp; width W (azimuthal 
extent of 3 time zones in the equatorial plane) v 15 R, and 


thickness b (the separation between the field lines from the 


poleward borders and equatorward borders of the auroral oval 


mapped to the equatorial plane) ~ 5 Re. Then the magnetic 


energy associated with the current flow in the volume gives 
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Figure 64 
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rise to the inductance Ly in Figure 62 and the kinetic Siege 
associated with the convecting plasma within the volume is 
represented by the capacitance Cu. The dissipative Pedersen 
currents flowing in the ionosphere are equivalent to thie 
dissipative current flowing through the Tonospheric resistance 
Ri. The high number densities of cn in the F layer and the 
changing electric field due to the transient response of the 
circuit gives rise to the displacement current through the 
capacitor Cy. Thus, the equivalent LCR circuit is based on 
the geometry of the current flow. 

Using the expressions developed by Rostoker and 
Bostrom [1976], we can determine the values for L, C and R 
and compute the resonant frequency of the circuit to appraise 
the viability of LCR circuit in representing the magneto- 
spheric circuit. The inductance Ly can be evaluated by 
estimating the magnetic energy associated with a current 
Uniformly distributed in the volume of length &.»% 15 Res 
with W. 15 Re and thickness: b ~ 5 Re. By equating this to 


the energy stored in an inductance carrying the same total 


eurrent,. ViZza; 9 
wbx 1 


2 
B ON eS oae 


we find that Ly TS <u? TH 


The capacitance represents the kinetic energy of 
the convecting plasma. Assuming a uniform plasma velocity 
Vere [7 <x 10° m/sec (y = p where E ~ 0.33 mv/m and B wv 20y 


in the equatorial plane) in the same volume as used in the 


Ym DFGonTR ony bas ‘$3. swe at 4) sonsiaubat 
2 omy Tay odd ning Pw smzerq patbaayeos aay sari 

en Pee ovis gteetb oat, eh) Sonsttosas> ciel ud b 
ent 0} sas Taviupe ate syadgzanot oils we ewnar 


S2nbs Ses atredgeonty ant gums Bhiwor FHSWIMD | 
odt bie 19¥BE 4 wad at *o: +o eatg tenet sedman Meta 
ant 0 s2nogea% insTengy? ont ot ‘aub bhatt. atisoele 2 
ght Hpuovds © PASTY. Jemsoetq2ib eat os s2iy zovig 

vo beesd et +tuaNts B04 sno fev hips and < eu 48 bia 
ee | wot trey odd 40 ‘ite 

bis. 1916S 08 ud boqotaveb, enotagorgxe gitd. nia . a 
1 bine 2) (gil ot asu! TAY ait animvesab oe ow Lavery m 
S70 g6 Pruantp ont 4b, oneupsat taanoeay edt. atngs 


-o snigatt ang anid aseeygen nf diworto A331 to. nee 


vd bospulays: od Nb? ce aonedoubat oat 5 hag 
J SY Hes a bess toozes \orsne- a idenesm and pats 


eget ef co a gpner. io emu fov orig nf basndagath, yim 


act eine ‘ontzevpa ya 


is 


rye o re apends tir bas itm we 
istos ome2 eat “ou senna ns nt bevode | mm 
“ae | Baa 


Sars 


Stag 84 | 
Roe Bu ap Ea ‘i #b oa Wt ‘a Po 
; \ , i | ! i | L } 


soe Ss - an sony batt 
xn ernst tosqss oat hn - 


io \pysng ottenia ‘ays ssnges | 
a amen tg” ontgasvnoa add 


Aad paves mit tay 8 


yOS' @ bas m\vm Sf. o f A Boo = y) o92\m toy x - Pay : 


ait af bee 26. “omy Fou anise oho at tase ratrogoupa: oe at 


152 
case of the inductance and, equating the kinetic energy in 


une "region to: the ‘energy stored timia capacitor’) ai zo, 
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= Dae =i) SENS WOE eae 4 
Weide pe an RPM Ry 
we find that Cy wee dO? Fr where we haveoused ou 8 Xx fo 
kg/m> Chor n 5 cm?) and ¥o © 10" vel tel. | 


For an F-region of thickness 100 km centered 
on ~ 300 km altitude, the ionosphere can (neglecting end 
effects) be represented by a parallel plate capacitor with 


capacitance 


where the dielectric constant of the plasma is 


ls 
n Mm, 


o being the mass density of the ions, n the number density, 

m. the molecular weight of 0° times amu and B the ambient 

magnetic field. Hence, « % 3 x 107° F/m for n wv 3 x 100 > Mee 

For A ~ 500 km x 3 time zones and d v 100 km, Cy nur OER 
Assuming an ionospheric rEsStTStivity OT 

i) aes is Qm averaged over the region of the ionosphere 

from ~ 110 km to 160 km, the ionospheric resistance below 


the F-region is estimated to be 
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From these values, it. can be seen that the. cir- 
cuit shown in Figure 62 can further be Simplified to a 
circuit shown in Figure 65. The values for the circuit 
_ parameters are summarized in Table 6. 

As mentioned earlier, the energy source is the 
kinetic energy of the plasma drifting through the dawn 
sector. The kinetic energy flux of drifting plasma through 
a surface in the dawn meridian of area of length 2 v 15 Re 
and width b wv 5 Re perpendicular to the flow is 
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mech | 


The electric power dissipated in the auroral electric 


region is 


Los omit aie IE,d 


Pore, @ 20 mv/m [Mozér and Lucht, 1974] and d ~™ 500° km and 
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ee 10° A, we obtained 
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The requirement of Peon > Pare yields V > 20 km/sec. 


Such velocities are consistent with the velocity in the 


equatorial plane near L ~ 10. Near bo 10;B ~~ 20y , 


bo 0.33 mv/m and Y = Fi uw 17 km/sec. «Thus, ‘the kinetic 


energy available in the plasma drifting past the dawn 


meridian is adequate to supply the energy requirements 
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TABLE 6 


Equivalent Circuit Parameters for Model Current System 


Inductance of generator region 


Capacitance of generator 


Capacitance of F-region 


Resistance of Ionosphere 
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associated with the ionospheric electrojet system. The 
equivalent circuit of Figure 65 is therefore valid. 

The model current system presented in Figure 65 
moeenrectively a steady state system. The results of our 
analysis of Pe 5 pulsations suggest that these micropulsa- 
mons are excited in conjunction with changes in the size and 
intensity of the electrojet system. We therefore wish to 
Study what occurs if we perturb the current in equivalent 
circuit shown in Figure 62. A sudden change of ¢ in the 
circuit will cause a sudden change in the current flow. This 
will set up an oscillation in the circuit whose frequency 


is that for the Li Cy DOwclOn— Ot mbike CT RCW IG ss Vilize 


For the equivalent circuit parameters shown in Table 6 this 
frequency is in the Pc 5 frequency range, Vil Zee AOS Lay any al 


and Cy # 3 x 10° F we obtain 


fue biz 


which is at the long period end of the Pc 5 spectrum. 
(Higher frequencies may be excited by exciting oscillations 
over a lesser latitude range or within a narrower electro- 
jet.) 

The oscillations we have been discussing above, 
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Figure 65 
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be imposed on the ionospheric portion of the circuit which 
consists of a pes istanke (R,) and capacitance (C,) in 
parallel with one another. Of course, the barbenteenbough 
Cy will lead the current through Ry by90° which provides the 
out-of-phase component of the perturbation magnetic field 
necessary to create a polarization ellipse for the pulsation 
in the horizontal plane. Were it not for the F-region 
displacement current, an oscillating current system could 
only lead to linear polarization for the pulsation. 

The model we have discussed above presents Pc 5 
aS an impulsive wave train, which belies the term “continuous" 
given to this class of pulsation. However, as we have 
pointed out earlier, the Pc 5's we have observed appear to 
be a series of impulsive wave trains which often blend 
together to give the impression that they are continuous 
in nature. The precise behaviour of the pulsations is 
controlled primarily by Cu and RT: Clearly, the more energy | 
there is in the generator region, the larger will be the 
effective capacitance of the equivalent circuit and thus 
the lower will be the lowest frequency in the Pc 5 spectrum. 
In addition the RC-time constant of the circuit will deter- 
mine whether or not the Pc 5 pulsations appear as a quasi- 
Continuous chain of oscillations or as a series of separated 
Again. it-is’clear*that the- more’ energy 


damped wave trains. 


there is in the generator region (i.e., the more energy there 
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is stored in the magnetosphere) the longer will be the 
meaeime constant for the circuit. This effect is reflected 


by the fact that while one might expect C,, to increase witn 


M 
increased energy storage in the magnetosphere, the conduc- 
Cavity 1s controlled primarily by solar: UV" in’ tHe morning 
Sector so that Ry will be relatively unchanged unless 
electron precipitation becomes extraordinarily high. The 
overall effect of the changes in L, C and R with changing 
levels of energy storage in the magnetosphere is that for 
lower energy levels the low frequency component of the Pc 5 
spectrum will be higher and the oscillations will appear as 
impulsive damped wavetrains while for higher energy levels 
the low frequency component of the Pc 5 spectrum will be 
lower and the wave trains will tend to merge with one 
another. 

From the above considerations, we would like to 
present the following model for the Pc 5 pulsations. 

A change in IMF (interplanetary magnetic field) 
and/or solar wind velocity produces changes in the convective 
velocity of the drifting plasma inside the magnetosphere. 

As a result, the energy available to drive the current 
system is changed. The adjustment in the size’ of the 
electrojet leads to spatial oscit lations and= the ciranger ih 
current intensity leads to intensity variations.” We spacrall 


and temporal oscillations of the electrojet combine to 
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generate the magnetic pulsations observed on the ground 
(see Figure 61). 

The actual currents responsible for the generation 
of Pc 5 micropulsations form a tithes editene nares system 
with currents flowing downward into the electrojet in the 
pre-noon quadrant, westward through the ionosphere as part 
of the westward electrojet (Hall current), and back into the 
magnetosphere at dawn (because of the conductivity discon- 
tinuity near dawn due to enhanced photo-ionization in the 
day time). The entire system is enclosed by the Birkeland 
current sheet pairs closed in the ionosphere through south- 
ward flowing Pedersen current (see Figure 5). Because of 
the relatively short longitudinal extent of the portion of 
the westward electrojet associated with the micropulsations, 
ground based magnetometers should be able to detect pertur- 
bations from both components of the three-dimensional current 
System. The contribution from both components of the current 
system is well illustrated by a series of latitude profiles 
shown in Figure 66. Figure 66 consists of six latitude 
profiles, each about 2 minutes apart. The borders of the 
convection westward electrojet are indicated by arrowheads 
at the bottom of each frame. It can be seen in Figure 66 
that AY steps of the latitude profiles on the left column 
wirem. profiles A, C and D) were located at lower latitudes 
than the AY steps of the latitude profiles on the right 


column (i.e., profiles B. D and F). The AY profiles of 
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latitude profiles B, D and F indicate the net downward 
field-aligned current flow across the electrojet while the 
meeproriles of ‘latitude profites ASC and E shows the net 
downward field-aligned current flow and the ‘end effect' 

of the westward electrojet portion of which is diverted up 
the field lines near dawn. Note also in Figure 66 the 
periodic spatial movement of the equatorward border of the 
electrojet. In order to gain a better insight to the 

changes produced by such a current system, we employed the 
techniques of Kisabeth [1972] to construct a model three- 
dimensional current system utilizing both the Hall current 
circuit and the direct current circuit involving the Birkeland 
current sheets. In the model a current system with 4° of 
latitudinal extent and 16° of tonecckianere extent was 
utilized where the Hall and Pedersen currents were distri- 
buted uniformly within the ionospheric electrojet region. A 
current of 10° A was employed and the total Hall current flow 
was set equal to the total Pedersen current flow for 
demonstration purposes. 

We show in Figure 67 the latitude profile taken 
along a meridian 4° to the west of the central meridian 
(which bisects the electrojet). The arrow heads on the 
Figure mark the borders of the electrojet. Note that the 
AY profile does not have a positive step going from south 


to north as is evident in the observed latitude profiles 


DES ay 


ae 
<a ees 
ha a) OP hea 


Mi \ eat Me 


“soothe bas! odd ‘be. inh “$n3} 


gil bodaey ib ‘ei stot dw to ‘nD EH & 5} tite [BY 
‘edt 38 yupia ni oats ston ; a sian eontt fai 

Sidi to Tsbred brenqatauns ond 40 Fngmsvon Fetsteqe at 
ary: ag Silp Fe re or “Wesied) 5 nitsg be Aebxb- ab , 
ads beyorgme SW agave daar 6 ‘Howe “ae weouhenytl 
“ast Fsbom 6 “touad nga, oF {svel} svoueety te) 291 
“ten ee ire ont dtod oats tt tw. modaye. $a! 7 
basladv(a edd. entyi Foun ftugt fo Fashtus saontiy oft 
ee - Avtw mas eye tiene 6 Tebom ous al 
a5W inate Ton tbud ignot 40 at bes. ere 
-fydetb: a96w eines 1u9. Heevebs bas TTeH odd otedw 
A notgss Jatargsasie, diradqzonot oi ntagiw <ensotien 
wot? FA SyANS, ‘Thee. fistod edd bas bevotqna: 2am, i Ponte a 
Tot wort grea, nge1sb9d Teves oat ot fsups 4 


UF, pte 
iY fei LA 


Mar 


ale mh ie | -2s2oqrug not 
noass gt re a As Fis af wote OW! 


a 


git ‘no ee 


¥ ‘ . 4 f 
and tent ato aBhordos Y [99 


ASU 2: ox? vaio pre vid feed is! evs ton 290b afi ya. " 
zat tong sbud ise! bevragay {fd f ettae? eT) = a - 
f f ; a) i! 1) 


162 


79 e4nbly 


(No) ‘LV1 SILJNDVWO39 
cc Pee Slee <°) SE eee eee: 


° 
®%@e50° 


1seeeeeesa' 


VN I'O= LNSYHND 340d WOILSYOSHL 


(4) SJGNLINSVW dI3ls 


: 


LIdnLe O\ 


Sof 


YL (oi) 


faa 
Set 


CEOWY CEL! 


| a hee ) j mia, 7 


| ‘on? (era warn 
9B 6 ada 
ee ay hina a dH. 


2 
) as y yy A } 


163 


Perause net downward field-aligned current flow was not 
accounted for in the theoretical profile. 

The electrojet borders were allowed to oscillate 
Spatially with an amplitude of 1°-of Jatitude.-and.the current 
was permitted to oscillate in strength with an amplitude 
equal to 33% of the total current flow. Figure 68 shows 
the differential profile for one quarter cycle of the 
Oscillation. The similarity of this profile with the 
observed differential profiles shown in Figure 69 is striking. 
Mme dataiused to construct the profiles in Figure 69 were 
4-100 mHz bandpass filtered to show the magnetic pertur- 
baicions of micropulsations: at ~ 5imHz. /The two profiles 
are about one quarter of a cycle apart. The polarity 
pattern can be seen to relate to the relative phase changes 
of each component across our line of stations as shown 
earlier in Figures 48 and 49 which indicate that the D- 
component normally does not experience a pronounced change 
in phase across the electrojet while the H-component ex- 
periences a 180° phase change across the selectrojet and. the 
Z-component suffers pronounced phase changes at the borders 
of the electrojet. An additional phase change in the H- 
Component is seen at v 74° in Figure. 68... Bins: 1s noted 
statistical feature of the data but is seen on some 
occasions (viz. the early hours of Day 16, 1972 shown in 


Figure 36). Of course, our choice of Hall current and 
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Pedersen current intensities is somewhat arbitrary (although 
the chosen values are quite reasonable). Variations in the 
current densities coupled with changes in the positions of 
the points of observation with respect to the current system 
will cause different relative phase patterns to appear. 
However, it can be seen that our model provides a aoe tone 
explanation of the relative phase changes of the components 
of the disturbance across the electrojet. 

The above considerations explain the observed 
relative phase changes across our line of stations. However, 
if there is no phase difference between the components at 
a given location, one would only observe linear polarization 
with the "lines" tilted in the first quadrant (components in 
phase) or the second quadrant (components 180° out of phase). 
However the ellipses ‘open’ up and switch in the sense of 
polarization across the stations, implying that there 
exist phase differences between components at a given 
location. We shall seek to explain these effects as follows. 

First. let us consider the polarization .in the 
H-D plane. The contribution to D comes from the direct 
"ohmic" current and a displacement current flowing in and to 
some extent above the F-region. The displacement current 
provides the out-of-phase component to D, which results in 
the total D perturbation being out of phase with H. The 


displacement current arises from a separation of the positive 
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defined as 
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where D is the displacement, E the electric field, e the 
dielectric constant, o the mass density of the plasma 


(= Nm. where N is the number density and m, = mean molecular 


weight X amu) and B the magnetic field. Using nominal 


values near the peak of the F-layer (i.e., mean molecular 


weight of 17, N~ 2.x 10° cm°, B~ 55,0007, and a conduc- 


: sm7!) and considering a period of »~ 200 


Bivity of" 2"x 10° 
sec (f ~ 5 mHz), we find that the displacement current is 


Mabout’ 1/10 "of the ohmic ‘current: 


5 ab 
fa) Eee 
2 


Thus, the displacement current has a small-but’significant 


contribution to the D perturbation. 


Since current going through a capacitance leads 


the current going through a resistor in a parallel RC 


circuit by 90°, one can express the perturbation in D 


observed on the ground as 


cos wt + Dy cos (@t + 5 (16) 


where Do includes the contributions to the magnetic pertur- 
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bation due to the ohmic component of the ionospheric current 
flow and the field-aligned currents (and is in phase with h) 
and Dp is the contribution to the magnetic perturbation due 

to the displacement current. Equation (16) can be expressed 
as 


D = D, cos (wtp to 8) 


where Do = D 


— DP 


6 =. tan 


6 is then the phase difference between the D component and 
the H component. It can be seen from Figures 68 and 69 that 
the H and D components vary in the Same sense south of the 
AX’ cross-over. ‘Thus’ 'D' leads H south of the AX’ cross-over 
producing counter-clockwise (Gey polarization in the H-D 
Orlane.’ North of’ the AX cross-over, H and D vary in the 
opposite sense. Hence, the phase difference due to the 
displacement current causes H to lead D in phase, producing 
the observed clockwise (CW) polarization (see Figure NA Ie 
The earth has a finite conductivity that causes 
the induced fields to be phase-shifted with respect to the 
inducing or source field. The phase shifting yor the: inter= 
nal field contribution results in complex Lorca ly 1 es 
inducing plus induced) field components at the surface. 
Nopper and Hermance [1974] have shown that the induction 


effect of the earth causes phase changes in H and Z with 


H leading Z. Again, from Figures 68 and 69; one can: see “that 
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as one moves from south to north, AZ and AX vary in the 
Same sense (i.e., H leads Z resulting in CW polarization 
in the H-Z plane), and then in the opposite sense in a 
rather limited latitudinal range (i.e., Z leads H resulting 
in CC polarization), and then in the same sense again (CW 
polarization) and finally in the opposite sense again (CC 
polarization). This agrees well with the observed senses 
of polarization in the H-Z plane (see Figures 37, 40 and 43) 
although over some intervals the CC regime in the middle of 
the CW Peat otis was not observed possibly because its limited 
extent made it difficult to detect using our stations. 

Atiasgiven station, Y, beads iH due se the ve ttec uso. 
the displacement current and H leads Z due to the effect of 
induction. It therefore follows that D leads Z. Again 
referring to Figures 68 and 69, moving from south to north, 
Pee can see that AZ and AY Vary in the ‘same sense {1-e.4 D 
leads Z resulting in CW polarization in the D-Z plane), and 
then in the opposite sense (i.e., Z leads D resulting in CC 
polarization) and finally in the same sense again (i.e., CW 
polarization). This agrees with the observed senses of 
polarization (see Figure ye 

The above considerations explain the observed sense 
of polarization in three orthogonal planes. The orientation 
of the major axis of the ellipse is due to the relative 


magnitudes of the components and to the size of the phase 
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difference between the components. A more detailed study 
on the polarization characteristics would undoubtedly 


yield additional information on the subject. of induction. 
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CHAPTER V 
'Pes4 CACTAMITY..UN TBHE-MORNING~SECTOR 


5.1 Observations 

As mentioned earlier in the first chapter, the 
occurrence of'Pc 4'giant pulsations is rare. We have only 
noted four events from almost five years of data recorded 
at our ions. They occurred on Vay 3475 1971, Day) 2o3;, 
254 and 255, 1974. The events all occurred in the morning 
sector. 

The magnetograms both unfiltered and filtered 

(with a 1-20 mHz digital bandpass filter) for the events 

are shown in Figures 70, 71 (Day 347, LOLA iee roe o. Coay 
PNG 74). 144075 (Day esd so 1974 ar omand 77 Chay 255) 1974). 

It can be seen from Figures 70 and 71 that the 
'Pc 4' pulsations for Day 347, 1971 are highly localized in 
latitude maximizing at MCMU and MENK and their highly 
sinusoidal wave forms are barely detectable at the northern 
stations Tike CAMB and CONT. “The-eventwof=Day.347 5. 1971 
lasted for about half an hour and appears as a single 'beat' 
or amplitude modulated signal. 

The events recorded in 1974 recur for successive 
days at intervals of about 94 hours; such a periodicity has 


been observed in~the past by Sucksdortff [1939]. The 
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unfiltered magnetograms for these events indicate that the 
magnetosphere was in a quiet state. The substorms that 
occurred around midnight were of small magnitude which 
hardly exceeded 150y. The filtered magnetograms show that 
the wave trains appear in 'beats' and that MCMU appeared 

to register the ‘Pc 4' wave train before the other stations. 
The 'beat' feature of the pulsation was most pronounced in 
the D-component at MCMU for the Day 254, 1974 event. The 
Pe 4' pulsation was obvious only at MCMU for the Day 255, 
1974 event. This event seems to occur in two separate 
bursts, the first one beginning at 0800 UT lasting for 
about half an hour and the second one beginning at 1000 UT 
lasting for.about an hour. “It Ts interesting to note, from 
“agqure /6, that the’ onset of substorm activity as seen in 
the midnight eector at the station of inuvik (71.9 W ‘geo- 
magnetic) coincides with the termination of the 'Pc 4' wave 
trains which were initiated at ~ 0805 UT. The substorm 
activity persisted until ~ 1000 UT at which time the 'Pc 4' 
activity reappeared. The Pi 2 activity) av TOW Tati tude 
stations were also evident during the interval. This 
behaviour of the micropulsations is extremely interesting in 
that it suggests that conditions of low activity must pre- 
vail in order that the source mechanism for” Pe a Rg arr 
It should be pointed out that 


pulsations can function. 


Rostoker and Falthammer [1967] noted that Pc 5 events tended 
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to occur in the decay phase of magnetic storm activity and 
that any Pc 5 activity in progress seemed to be terminated 
at the onset of a magnetic storm. Thus, there is an 
apparent /vendency for both Pc 5 and Pe 4. giant pulsations 
to be associated with decaying or low level activity in the 
magnetosphere, while the source mechanisms appear to be less 
effective during periods of increasing magnetospheric 
activity. | 
To study the frequency content of the 'Pc 4', raw 

power spectra were computed and are presented in the next 
few diagrams. 

| Figure 78 shows the raw power spectra for Day 347, 
O71 in ithe “intenval of 1745 and 1845) Utyy .The erratic 
appearance of the spectrum is due to the fact that the 
spectra were unsmoothed (see Figure B2 in Appendix B for the 
appearance of the peaks resulting from smoothing). The 
frequency of 'Pc. 4" is revealed by the spectra to De in the 
10-12 mHz band. The peaks in the 10-12 mz band seem to 
shift from high frequency to low frequency going from south 
tomnorchs ‘The double peak. structure 7s evident at) SMI, 
FTCH, MCMU and MENK. The double peak feature might be due 
to statistical fluctuations of the estimates (see again 
Figure B2 in Appendix B for the variations of the confidence 
interval and the bandwidth) and might not be real. However, 


if it is real, then the ‘beating’ appearance of the ‘Pc 4! 
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pulsations on the magnetograms can readily be explained as 
due to the addition of two waves with nearly the same 
frequency. The northern stations such as CONT show large 
powers in the low frequency band of the spectrum and some 
high frequency peaks across the Pc 4 band. This is 
associated with the convection westward electrojet noted 
at the northern stations of CAMB and CONT for this time 
interval@as shown. in Figure 79. The ‘Pe 4° peak in the 
spectra which is evident at the southern Stations well 
away from the electrojet is evidently not related to the 
convection westward electrojet which has been shown in the 
previous chapter to be the source of Pc 5 pulsations. 

The spectra for Day 253, 1974 in the interval 
between 1330 and 1430 shown in Figure 80 indicates that 
the 'Pc 4' pulsations have frequency in theakl.5=13:55 MHz 
band. The double peak features are evident at all the 
stations except MCMU which shows multiple peak structure. 
The most dominant peak is centered at 12-12.5 mHz. 

Figure 81 shows the spectra for Day 254, 1974 in 
the interval between 0945 UT and 1045 UT. The ‘Pe 4! 
peaks are clustered in the 11-12.5 mHz band. The double 
peak feature is evident at most stations, cler vl... Mic. pedn 
being the most dominant one. The most remarkable feature 


of these spectra is the complete absence of powers in any 


other frequency except the ‘Pc 4’ frequency in the D compo- 
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nent at URAN and MCMU. This implies. that the signal is 
virtually noise-free reflecting the quiet state of the 
magnetosphere and that the signal is highly monochromatic. 

; The spectra for Day 255, 1974 in the invervea | 
between 1000 UT and 1100 UT are shown in Figure 82. The 
"Pe 4) peak centered At oF MZ aS OY pronounced at 
MCMU. The double peak feature is evident in the D-component 
at MCMU. Note also the complete absence of powers in any 
other frequency except the 'Pc 4' band in D at MCMU. 

To study’ the characté@ristics ofthe 'Pc 4’ 
pulsations, polarizations and phases were determined using 
the smoothed Spectral estimates. The powers in three 
components, the relative phases across the stations and the 
polarization ellipses were plotted as a function of 
invariant latitude or dipole longitude for the 'Pc 4' spectral 
component for each event. The results are shown in the next 
few diagrams. 

Figure 83 is the latitude profile for the pulsation 
characteristics atel1l.2 mHz for Day 347, 1071 in the inter- 
val between 1745 and 1845. The latitude profile of the 
power findicdtesthat the evant-is, highly localized in 
latitude. The H power profile suggests that the peak power 
probably lies between MCMU and MENK. The power drops off 
rapidly to noise level in about 4° away from the peak. The 


relative phase plot shows that there is rapid phase change 


of Tenote ans tons 29 tfgat’ eat now ee aaau 
| add to ata z getup ans pntazetien ‘sent-seton y 
oH smo soon vider ef vente ont teat BAS | 


Pewrsani ons mi AVOE .des. yd ‘et BtsI9q2 tye i 
0Ott bas Po ngaM 


BAT +188 eet at avwode ats. 1g 0 
ve. beonwonorg xine af xh err bevednaa 
Swenaqaos- <0 add af drektve ef erwt5o? 4699 stduop om 
| es. i 21swoq To sonbeds dof qmos en? dels aton Be. ; 2 ik hi 
UMOM te 0 at bosed “8 34° x tq99x9 conoupent 

‘S 39" Shy We rea idetvatosieda eds ybuta ot a ai 

paren, bantmrad 9b 5194 292a8fq bas anotisstystoq vem ey iy 
eerdd af evswoq oAT .potemhg2s: tantooqe be | 
‘suid bos enotist2 aft eeoras: eseeny eviselay ONT. om qi 
to notssanut 6 26 beanerg 81s eseqtt Ts Lorde: 
, fsngoaae ‘DS 99' gat a0T shut remol afogth 10. obustael 
| Kon ‘ia nf ware sv6 astueew sat Saave ase, tad aa 
cameo ond 1a of rot suing FES 2t 8. pre ‘ ae. ¥, 
cesar ads nt FVer . WHE ysd 10% stm $. ‘7 46 eottety nF 
4 odd to of toxg obug ttet eat: “.anst brs anti ne 
id digit ii besttssor xe tM ab ; . 9 ont tedt astestoad 19 
awog A890 i tend eseavue at itor 49K0q Hq eat 7 at 
ee A vary v9Mod onT wan bas uMoM dona estt Psdon 
dbeq sit mont yews 8 Sunde at fevat szton oo yrbtasy 


y 


Bet g2stq bigs et syede todd aworte tofq ses ovitstoy 


; : 
“ iis 
Th Sey 
es 
« 
oy 
iT 
aes 7 
3 
4 a 7 


— i 


zg eunbl4 


7HW NI O34 - tHW NI 0344 zHY NI 03u4 
QISIVIGIZITIONG @LHSHE STO SIIPIEIZIIIOLE SL OSH ES LO SISIHIEIZIIIOIG6 BLUSH EZIO 


fe H 
(LOOJb-OO0TL £61. SS2 AUG 


0° ot 


THY /oL 


a 


-_ 


eg eunbLy 


Sj WW1 


LVINI 
ZHW Z'Il (ShBl-Sb2I) IZ6I “LbE AVA 


Oh a 


Nn 


Nfeaweseccccceee ecccoss 
x 


(zHwWA4) YSMOd 


| ; ohi! ace 


190 


between MCMU and MENK in H and between FTCH and SMIT in D. 
The phases aut higher latitude stations were not plotted 
because of their low powers. The polarization ellipses 
projected in the H-D plane indicate that there is a switch 
in the sense of polarization between MCMU and MENK from 
Uwe sense dt polarization.to..6€’sense of polarization. In 
ocing trom SMIT to .LEDU, the ellipses if H-D plane are 
tilted towards, the first quadrant except at MCMU where 
it is tilted towards the second quadrant and inclined 
towards the D direction indicating the Vamge power in D. 
In going from SMIT to LEDU, the ellipses in the H-Z plane 
are aligned predominantly in the H direction with the 
kcceeae of MCMU where the ellipse is tilted towards Z. 
The sense of polarization switches from CW to CC 
between MCMU and MENK. The sense of polarization in the 
H-Z plane for the stations where the effect fomeirc = 
pulsations is felt is the same as the sense of polarization 
in the H-D plane. This behaviour is contrast to that of 
De & discuscetoim.the last chapter. Forgthe Pc 5 pulsations 
whose origin lies in the electrojet, the sense of polariza- 
tion in the H-Z plane is opposite to the sense of polariza- 
tion in the H-D plane. 

| Figures 84, 85 and 86 show longitude profiles of 
the characteristics of the'Pc 4'giant pulsations recorded 


on Day 253, 254 and 255, 1974. The data from MCMU fare 
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DAY 253,74 (1330-1430) at 12.5 mHz 
DIPOLE LONG. . 
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Figure 84 
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DAY 254,74 (0945— 1045) at ll. 7 mHz 


DIPOLE LONG. 
HAYR SMIT MCMU .URAN 


- Figure 85 


= a our i y | ‘ 


ot 


a era ag . 
. alindhgasee vehi a YY 


\ aU, ‘ tt Nay Ue ME ‘ nS 4h \ ee 
‘oo ee Ooi warn yj ai j mo ot } zy 7 mii, hs, = 4 7 a i 
| ay ee ta) a (Nah aaa. 1h ‘ ee a ae i make ie fan 
i ; ai >, ’ ; 5, 2 wag we 
Va 


tit 
oN 


Pe a 
Biel anne be Pr 
ma! te : hes aad 


i 


DAY 255,74 (1000~1100) at 1.7 mHz 


DIPOLE LONG. 
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Figure 86 
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included on these profiles to allow comparisons to be made, 
although it must be remembered that MCMU is ~ 3° equatorward 
of the other three stations and is nearly due south of URAN. 

On Day 253 we see the power at 12.5 mHz over the 
Universal Time interval 1330-1430. The power is idgest 
at MCMU, but tends to vary little in longitude over the 
~ 15° span of the coverage. The change in sense of polariza- 
tion in the H-D and H-Z planes between MCMU and URAN is 
consistent with the source region being centered between 
Preseetwo' Stations. lt as interesting to note thatethe 
senses of polarization in the H-D, H-Z and D-Z planes do not 
change across the longitudinal span of the array. However, 
there are outstanding phase reversals between URAN and SMIT 
and between SMIT and HAYR in all three components. We 
shall discuss the significance of these phase changes later 
if this "chapter. 

On Day 254 using the power at 11.7 mHz over the 
interval 0945-1045 UT we see behaviour of the micropulsation 
characteristics much the same as on Day 253, although the 
power levels are slightly lower. Again the peak power is 
at MCMU and there are phase changes in the H-D and H-Z planes 
between MCMU and URAN. Aside from a tendency for linear 
-poldr#Zatiron rm the D-Z plane at SMIT, the polarizations 


are unchanged over the longitude span of the array. Again 


the outstanding changes in relative phase in all three compo- 
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nents between URAN and SMIT are present but between SMIT 
and HAYR there is little relative phase change. 

On, Day .255,,using: the power at dberethz Oyvyerethe 
interval 1000-1100 UT, we again see behaviour of the 
micropulsation characteristics eee ke inOSe, for 
bus 254 and 255. or. this.day.thepsignaleis highly 
localized in latitude, with MCMU being closest to the 
center of the source region. However, the power levels are 
considerably reduced from those on the previous two days. 
The relative phase changes in longitude precisely parallel 
enose for Day 253. 

For purposes of comparison with previous workers 
we draw attention to the detailed studies of Pc 4 giant 
pulsations which have been carried out by Eleman [1966] 
using data from eleven high latitude stations in Scandinavia. 
He claimed that these pulsations were distinguished by the 
fact that the vertical component of the perturbation is 
almost as large as the horizontal component near the center 
of the disturbed region. In agreement with Eleman we find 
that the vertical component may have magnitudes comparable 
to that of H near the source region although the power in 
the D-component usually exceeds that in the H-component. 
Although he found for most of the events that the pulsations 
were polarized in the CC sense in the horizontal plane, 


on some occasions he observed the polarization in) te 
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horizontal plane to be polarized CW at the northern 
stations ang CC Bethe Southern stations. The latitudinal 
profile of the power and ellipses of the event we studied 
Suggests that the stations north of the power peak had CW 
ellipses while eediabns to the south had CC ellipses in the 
horizontal plane. He also contended that the region of | 
disturbance tended to remain fixed geographically during 
the lifetime of the ‘Pc 4' pulsation event. He was unable, 
however, to settle on any theoretical explanation of the 


source’ of ~Pc’ 4° pulsations. 


5.2 Generation Mechanism for 'Pc 4' 

Veldkamp [1960], after studying.a "Pc 4" event 
which had maximum amplitude at ~ 57°N geomagnetic, contended 
that the.observed period of 105.sec fits well with eigen- 
oscillations of period ~ 100 sec predicted by Obayashi and 
Jacobs [1958] for field lines penetrating 57° geomagnetic 
latitude. Annextad and Wilson [1968], after studying ‘Pc 4' 
events recorded at College and its conjugate, Macquarie 
Island, suggested that the 'Pc 4' is a result of the first 
even mode of transverse hydromagnetic wave propagating along 
field lines to conjugate points. 

it would be of interest to estimate the eigenperiod 


of the toroidal and poloidal oscillations for field lines 


emanating from L=5.4 which is the L value for MCMU where the 


amplitude of 'Pc 4' was maximized. We estimate the eigen- 


period using the results of Orr and Matthew [1971]. 
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Starting with the wave equation (4) for hydro- 
magnetic waves in a cold collisionless magnetized plasma 
derived earlier in Chapter II and assuming hice the magneto- 
spheric plasma is composed of protons and electrons with 


ie m 


a number density distribution n = neta) 


» where n and ie 
are the proton number densities at r and Fp respectively 
A Ge being the equatorial geocentric distance), Orr and 
Matthew derived an equation for axisymmetric toroidal 
oscillations and another equation for asymmetric guided 
poloidal oscillations. Both equations were integrated 
numerically for different plasma distribution defined by 
m= 2,3 and 4 and using the boundary condition that the 
alectinc -field,is szeresat;theseamth’s xssurfaces,) kheaperiiod 
for ae ieala field line with maximum equatorial 
geocentric distance My a L is given by 
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kc = wen the eigenvalue for axisymmetric toroidal 
Or 


oscillations or asymmetric poloidal oscillations 
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Using equation (17) for L=5.4 and the number 
densities of plasma in the plasma trough and inside the 
plasmapause in the equatorial plane (see Chappell, 1974), 
we obtain the following eigenperiods compiled in Table 7. 

| the eigenperiod T is in seconds and the frequency 
fase ttnumthiz.& lt can be: seen from the Table 7 that the first 
harmonic asymmetric guided poloidal oscillations for number 
deRsiities of 5 on 6 cm > and 10 cm” have frequencies 
ranging from 8.9 mHz to 13.6 mHz depending ore number 
density and the power law of density distribution. The 
observed frequencies of our 'Pc 4' events lie in that range. 
The first harmonic axisymmetric toroidal oscillations for 
ve 10 cm”? and the second harmonic axisymmetric toroidal 
Peiipaianns For Deange 100 cm also have frequencies very 
close to the observed frequencies. Equation (4) is a very 
simplified representation of the much more complicated 
conditions in the magnetosphere. As mentioned earlier in 
Chapter II, toroidal and poloidal oscillations are generally 
coupled together. The shear mode is decoupled for resonance 
condition. The above calculation is to obtain a rough estimate 
of the range of frequencies within which the experimental 
observation should lie. The results are thats the 
computed resonant frequencies seem to agree with the 
observed frequencies, thus giving some support to the 
contention that field line resonance might bepthe cause of 


the 'Pc 4' pulsations as it is definitely not related to the 
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auroral westward electrojet which is the primary source for 
Pe S*pulsations?°°Furthermores the reversal of sense of 
polarization in latitudes across the intensity peak is 
predicted by the field line resonance model of Chen and 
Hasegawa [1974a] and Southwood [1974]. The monochromatic 
character of the 'Pc 4' events as demonstrated by the power 
spectra (see for example, the power spectra in D for MCMU 
and URAN in Figure 81) suggests that the driving force is 
composed of nearly monochromatic waves. This possibility 
was solved by Chen and Hasegawa [1974a] and Hasegawa and 
Chen [1974] as a steady state case. The solution exhibits 
Surface wave away from the resonant field line and shear 
Alfvén waves near the resonant field line. The ‘Pc 4' events 
presented here therefore strongly suggest the excitation of 
shear Alfvén waves by a monochromatic source. The tilt of 
the ellipses in the horizontal plane does not agree with 
that predicted by Chen and Hasegawa [1974] and Lanzerotti 
et al. [1974] who showed that the major axis direction does 
not switch across the resonant region and lies in the 2nd 
quadrant. Our observations indicate that there is some 
switch in the tilt of the ellipses across the resonant region. 
This might be due to the possible effects of the ionosphere 
on the rotation of the ellipse orientations [Hughes, 1974]. 
Finally, we would like to discuss several charac- 


teristics of this type of pulsations which are unique. An 
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examination of the standard magnetograms and an estimation 
of the wave number m using the relative phase change in 
longitude (m ~ 36) indicate that the phenomenon is highly 
iMealtzed in longitude. “Pe 4. pulsations might bea 
manifestation of the quiet state of the magnetosphere as 
exemplified by the event of Day 255, 1974 where the onset of 
Substorm stops 'Pc 4' which reappears when the substorm 
activity terminates. 

The daily recurrence tendency of the 1974 events 
suggests that the pulsations might be related to some 
structure within the magnetosphere (such as the plasmapause) 
which is maintained over a significant period of time. However, 
the plasmapause extends over all hours in local time, and 
hence one might ask why the oscillation region is SO con- 
fined spatially. This leads one to look for some significant 
magnetospheric parameter which may suffer some transition 
in the local magnetic time regime occupied by the pulsations. 
It is intéresting to note that the electric field observa- 
tions of Mozer and Lucht [1974] (see Figure 63) indicate a 
re eae from westward electric field to eastward electric 
field in the early morning hours and near noon. Such a 
transition would be associated with magnetic flux tubes 
drifting equatorward in the region of westward electric field 
and poleward in the region of eastward electric field. Thus 


field lines mapping to the plasmapause may suffer shear flows 
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av the region of,electric field transitions and/thisamay 
well prove to be a region which can be perturbed should the 


Character of the convective flow be subjected to any change. 
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CHAPTER Vi 
SUMMARY OF RESULTS AND CONCLUSIONS 


In this thesis we have studied in detail the 
characteristics of Pc 5 micropulsation activity in the 
morning sector and we have documented several cases of the 
rare 'Pc 4' giant pulsations. The following summarizes the 


major findings of this thesis. 


6.1 Pc 5 Micropulsations 

One significant finding of this thesis is that 
there is no clear latitudinal dependence of the principal 
spectral components of the Pc 5 pulsations. The observation 
that altethe eee components peak within the same 
latitudinal range occupied by the convection westward 
electrojet rules out the possibility that Pc 5 pulsations 
are generated by the resonant oscillations of field lines 
for a reasonably normal density distribution of thermal plasma 
in the magnetosphere. This finding and the results of the 
analysis of the intensity contours and polarization parameters 
of the major spectral bands using the convection electrojet as 
a frame of reference point strongly to the three-dimensional 
current systems as a primary source for the Pc 5 micropulsa- 
tions. In fact, we have shown that the spatial oscillations of 


the borders of the electrojet occur at frequencies which are 


clearly identifiable in the Pc 5 micropulsations. 
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We have emphasized that Pc 5 wave trains appear to 
be composed of a series of impulsive wave trains which merge 
together to present a picture of continuous activity. The 
response of the Pc 5 micropulsations to rapid changes in the 
configuration of the electrojet lends weight to the conten- 
thon that Po 5 pulsations .should .be. thought of .as,.a..Sseries .o,f 
impulsive wave trains, each wave train being associated with 
a rapid change in magnetospheric parameters. In particular, 
the rapid poleward motion of the equatorward border of the 
electrojet are associated with marked enhancement of power 
imethesPec5dspectrums .iThe Pe.S oscillations therefore 
appear to represent the transient response of the three 
dimensional magnetosphere-ionosphere current system associated 
with large scale convection within the magnetosphere. The 
fact that the reconfiguration of the electrojet appears to 
occur in response to changes in the north-south component of 
the interplanetary magnetic field lends weight to our 
suggestion that sudden changes in the level of convection 
in the magnetosphere are responsible for the initiation of 
ringing of the current system (as proposed by Bostrém [1972] 
among others). 

Based on our observational results, we have 
presented an equivalent circuit representing the magneto- 
sphere-ionosphere current system. We have postulated that 


the current system is driven through the extraction of 
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kinetic energy from the plasma convecting sunward past the 
flanks of the magnetosphere, in manner discussed by Rostoker 
and Bostrém [1976]. A sudden change in the convective 
flow results in a sudden change of the convection electric 
field (which also drives the westward electrojet) resulting 
in a perturbation of the electric circuit. The perturbation 
will cause ringing at the frequency of the LC oscillations 
of the equivalent circuit. We have demonstrated that these 
oscillations would be in the Pc 5 frequency range for a 
circuit involving downward current flow near noon linked by 
the westward electrojet to upward flow at the dawn meridian. 
Finally we have shown how a combination of spatial 
and temporal variations of the electrojet can give a 
plausible explanation of the relative phase change across 
Our stations. Furthermore, by introducing the displacement 
current in the F-region and induction in the earth, we can 
explain the observed senses of polarization in the three 


orthogonal planes. 


6.2 °'Pe 4" Micropulsations 

In contrast to the Pc 5 micropulsations, we have 
found évidenée that thestPe 4' wdiant ‘pulsations tare snor 
related to the westward electrojet although they appear to 
be influenced by the initiation of substorm activity. The 
latitude profile of the power in the principal spectral 


component reveals that this type of pulsation is highly 
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localized in latitude and the large azimuthal wave number 
deduced from the phase difference in longitude indicates 
that the activity is also highly localized in veintghieuticn 
This feature might be related to the transition of the 
azimuthal electric field. The observed frequency of 'Pc 4' 
pulsations compares favourably with the eigenfrequency of 
the field lines penetrating the region of maximum amplitude 
for the asymmetric poloidal oscillations. The highly 
monochromatic nature of the event and the reversal in the 
sense of polarization across the region of peak intensity 
give support to the popular field line resonance model by 


Chen and Hasegawa [1974] and Southwood [1974]. 


6.3 Conclusions 

There are undoubtedly different source mechanisms 
for different classes of ultra low frequency micropulsations. 
There might even be competing sources within the same class. 
From the detailed analysis of the Pc 5 micropulsation events, 
we contend that the three dimensional current system asso- 
ciated with magnetospheric convection dominates other sources 
in its contribution to the low frequency band of the micro- 
“pulsation spectrum. However, we would expect that away from 
the influence of the electrojet during quiet conditions the 
field lines might be excited to resonate at their natural 
frequencies under certain conditions. The rare ‘Pc 4' giant 


pulsations may well be due to fieid lihe-resonances... ineiivys 
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it should be pointed out that long period micropulsations, 
in addition to being a monitor of the solar-terrestrial 
interaction, might provide a useful tool for the study of 


induction in the. earth. 
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APPENDIX A 
INSTRUMENTATION 


The equipment used to measure and record the 
fluctuations of the earth's magnetic field is basically a 
three component fluxgate magnetometer, whose output (in the 
form of DC voltage levels) is recorded on magnetic tape 
using an analog-to-digital recording system. The recording 
of data in digital form enables maximum efficiency for data 
handling through the ability to directly process the data 
using a computer. The basic components of the recording 
system are shown in Figure Al. Each system includes a 
three-component fluxgate magnetometer, analog-to-digital 
converter, digital tape recorder, WWVB time code receiver 
and the necessary external power supplies and regulators. 
Data were recorded on seven-track magnetic tapes at a rate 
of 1.56 samples/second (corresponding to a sample interval 
of 1.92 second/component). The records were timed by 
recording a WWVB time code signal for 2 minutes directly on 
tape at ‘7 5-hours intervals. 

All magnetometers, except the one at Leduc (where 
a Sharpe MFO-3 was used), were three-component fluxgate 
magnetometers of the type designed by Dr. P. H. Serson of 


the Geomagnetic Division of the Department of Energy, Mines 
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and Resources and manufactured by. £.D.A. Electronics Utd. 
The theory behind the fluxgate system has been discussed by 
Serson [1957] and the design specifications and results of 
tests performed on the magnetometer under normal field condi- 
tions have been presented in detail by Trigg et al. [1970]. 
The Sharpe MFO-3 magnetometer has similar specifications | 
except for better long term stability (+ ly/24 hrs at con- 
Stant temperature), a smaller temperature coefficient 
ian OU py/°C). and a broader frequency response (Dis Ge gUOM to De siize, 
-3db). | 

The sensing head of the magnetometer was mounted on 
an aluminum pole set in concrete at the bottom of a plastic 
container. The head was located where it was not likely to 
be disturbed by moving magnetic objects. The container was 
buried in the ground to reduce diurnal temperature varia- 
tions, the top being ~ 6 to 10 inches below the surface. The 
outputs from the X, Y and Z amplifiers are linear over the 
dynamic range (+ 10V) as shown in Figure A2. The phase and 
gain curves are essentially flat from DC to 0.1 Hz as shown 
in Figure A3. The sensitivity, resolution and dynamic range 
of the magnetometer are | Volt/WO0y. flys and: +. h000y 


(from baseline) respectively. 


The major problem with the fluxgate magnetometer is 
its sensitivity to temperature variation. The value of the 


temperature coefficient given by the manufacturer of the 
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Figure A2 
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magnetometer is less than ly/°C. However, the stability and 
drift tests completed by Trigg et al. [1970] showed that the 
fluctuations of the Y baseline (up to 20y) were much larger 
than expected and could not be explained by temperature 
drifts. Also a temperature coefficient of 10y/°C was 
| necessary to explain the variation of the Z-component 
baseline. 

The analog-to-digital system is based on the Redcor 
720 A-D converter, with multiplexer bypass line for WWVB, 
WWVB-to-data relay, and a variable stepping rate system 
coupled to a crystal “oscillator as shown in..Figune A4. The 
multiplexer is effectively three switches, one for each of 
the three Bamanwant channels. controlled by the sequencer. 
The sequencer controls the sampling order CAD ZED Ze eos. 
initiation of sampling of the multiplexer and the conversion 
of input voltage in the analog-to-digital converter. Before 
the signal from the multiplexer is fed to the analog-to- 
digital converter it is buffered. The input of the buffer 
can be coupled to the output of the multiplexer or the output 
of the WWVB receiver. This selection is controlled by a 400 
Hz crystal oscillator (with a maximum deviation from nominal 
frequency of + 0.005 per cent) and the associated divider 
circuit. The stepping rate was changed from 1.56/sec (sample 
anterval of 1.92 sec/component) to 50.0/sec while the WWVB 
time signal was being recorded. The relay switches from 


data to WWVB time signal on every eighth block. 
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The 12-bit A-D converter (the Teast significant 


bit was used as a flag to identify errors in the sampling 
order) output was formatted into 2 byte samples (each byte 
being 6 bits long) and recorded on a 7-track Peripheral 
Equipment incremental write tape deck. Each block contained 
6144 samples and represented 65 min 32.2 sec of magnetic 

data or 122.8 sec of WWVB time signal. Since the dynamic 
range was + 1000y and there were 1? digital bi tseseme 

Gamez ing anterval (amplitude sample interval) is 0.976y 
(20007/2!!) which is approximately the same as the resolution 
of the magnetometer (+ ly). 

The component outputs from the magnetometer were 
passed through two stage low-pass filters to remove frequency 
components above the Nyquist frequency (0.26 Hz) <toxavoud 
aliasing. The filters also served as buffers to isolate the 
multiplexer from the magnetometer by limiting the voltage 
input to less than + 12 volts. 

Because of the critical sequence of events that 
must take place to properly arm the unit for recording, an 
automatic restart unit was used to turn the system back on 
after a power failure. The auto restart has built-in logic 
to initiate control pulses to the logic board and tape 
transport in proper sequence to ensure the unit starts up 


properly. The WWVB time signal was also recorded immediately 


following the power failure to increase the. timing relvabi lative 
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APPENDIX B 


DETERMINATION OF POLARIZATION PARAMETERS AND 
PHASES FROM 3-DIMENSIONAL VECTOR TIME SERIES 


Previous methods for estimating the characteristics 
of micropulsations such as mean-periods, mean-amplitudes Ai 
polarization hodograms directly from the raw data Cite 
amplitude-time recordings of the perturbations) are often 
cumbersome and restrictive and the results are sometimes 
misleading [e.g. Pope, 1964; Egeland, 1965]. Recently, 
researchers in the field of micropulsations have adopted the 
cross-spectral formalism to determine the polarization 
characteristics of micropulsations. The approach is based 
on the quasi-monochromatic wave theory. .of phys uealmopeiesmas 
outlined in Born and Wolf [1959]. This technique has been 
applied by Fowler et al. [1967]. Rankin and Kurtz [1970], 
Samson et al. [1971] and Fukunishi et al. [1975] among others. 

From an analysis presented by Born and Wolf [1959], 
the polarization parameters for a quasi-monochromatic wave 
field at a given field point can be obtained in terms of the 


elements of the following 3-dimensional covariance matrix 


‘é eon ee 
x Ht Ve tee 
* ia 
ii ys 7 ue ae Wy 
Vaan mt ie i, F wa, . ay 
ith 8 f ai i is ’ Cv a fi 7 Lond a a 
poy veh dian 7 hae 


| sw ABT INAS HOTTAS seve a A | : i 
, ii <6 MORT 2320N4 


“9. of 6isb wey edd mort etioaiie este not : ee 

nedto ais (ehcisadiydi1aq ant to ‘Spb indied 4 utd -sbud ER 16 

: “pemitomoe evs eiTwean old bas ovitatagesy bis ‘om ne | 

¢¥finsooh [der (bel eg’ gbael .3qo% .e. oT e nied: 

ot } ber qobs even enottseluqors ti to brat oft nt ane se 28 

| hotyesivstoq aay onimretob od mei temvot lavtosqe-deo 

beesd ef dabovqqe sat eho tied hugess ta” to ane 6 ) 

26 2otiqo faoraytg to vroeny evew ‘Sid emoddsénod:tesup' 8 ay no 
need end sip taaser eit -feaery “Holl bas nves nt be ): 7 

/LOXOT] siawt das nilaea 18800) : $5" 18fwod co bat 

-2tondo enone hevery Te ‘ya. faatnbs vt bas Crveny Te 29 mozm 

| feeer] low bes avoe a dodiaaoia stators ns mova ote 

avew >tyenorrtoonom- hel up s 40 ebtonsteg notdes. oe) it gq odd. 

ond 70 aire ni bontsddo ‘od re ‘shiog biott novie 5 38 Wied 

xtaea son #15¥09\ Tenotensmrb=t pnt onde and Yo ethanol 


a 


j we ave ie Tick, boat a t ‘446m ae eho | og net tru qain ya al | 
Nin ARS a Bits au iat? | “eh anant ' q 
SM pind 'F:) HGR Maowe.g sigh Reine ae 


a] : ae 
? i ' i ; ; 
an : 
, ; 


ess F Al 


<B(t)B,(t)> BR CEI 8 LE) <B(t)B(t)> 
J=  |[<Bi(t)B(t)> <BL(t)BI(t)> <B.(t)B.(t)> 
<B,(t)By(t)> <B,(t)BY(t)> <B,(t)B_(t)> 


where the B's are the orthogonal magnetic field components, 
the brackets indicate a time average and the asterisks denote 
the complex conjugate. In general, the time-averaged 


elements can be expressed as 


" agen can: " 
Suna n= f B,(t)B,(t)dt 


for each Ja, element. This integral is evaluated, in prac- 
tice, by using the relationships among the covariance 


functions and the associated cross-spectral functions. This 


is expressed as 


’ te oleae Cade 
,5(7) 3 t! ; B.(t tT) j 
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where ee «5 the covariance function, t is the time lag, and 


Af is the effective bandwidth of the signals -When=t =a, 
55 is the auto power spectral density. When i # j, o,, 18 


the cross power spectral density. When t = 0, we have 
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Hence, the polarization parameters can be formulated in terms 
of the auto power and cross power spectral estimates. 

| The property that the real symmetric part of the 
cross spectral matrix is diagonal in the principal axis 
system is utilized to determine the apparent plane of 
polarization. The real part of the cross spectral matrix 


is diagonalized as follows: 


J4q 0 ) 

# |= ; 
R[Re JJR = 0) Joo 0 
0 J 33 


where R is the rotation matrix and we have chosen 017795075 33° 


The transform of the total matrix is then 
ai de eon 
Rok © = -idd5 52 1J53 


Gid5 a) oly geo. 


If the eigenvector associated with the minimum eigenvalue 


Cire! § J33) is used as the wave normal direction and the 


full matrix is rotated so that one of the new orthogonal 


axes corresponds to the wave normal direction assuming that 


the non-plane polarized vector series 1s Wyandom, tne 
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cross spectral matrix of the plane polarized signal is 


given by 
(5, a J33) iJ45 0 
-idi5 (J55 J33) 0 (Al) 
0 0 0 


Hence, the degree of polarization can be expressed as the 
ratio of the apparent plane polarized signal intensity to 


the total signal intensity as follows: 


where we have used the factathat.the trace of the matrix is 


equal to the intensity. 


At thasepome, cle ?-dimensional analysis of 
Fowler et al. [1967] and Rankin and Kurtz [1970] can be 
applied to the matrix representing the plane polarized 


signal. The plane polarized signal is composed of a vector 


moving randomly in the plane jn which the signal is polarized 


and another vector whose two orthogonal components are 


completely correlated. Hence, the matrix (Al) can be 


expanded in the form 
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where P is the matrix of the completely correlated vector 
series with det[P] = 0 and U is the matrix of the random 
plane polarized component. Solving for U (the eigenvalue 


or characteristic root of the matrix J), we get 


—" 


7 9 ; t/zZ 
BPI TT f SLC, ie ceina) oe eet 


QB 
The root with the positive sign is rejected since it gives 
negative values for the auto-powers of the polarized 
signal Pay and Poo: 

Hence, the matrix of the completely correlated 
vector series can be expressed in terms of the original 
| plane polarized matrix J' and the ratio of the completely 


coherent signal intensity to the plane polarized signal 


intensity given by the expression 


(J5 4-933 


(d1,-J54) * (9-933) 


=U) “F “(Jt ,-d44-U) 
“ee 22 toe 


can be determined. 


To interpret the matrix P in terms of the quantita- 
tive polarization parameters, consider the following 


mathematical illustration 
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Expanding the above equation gives 
ye 2 Re P 2 
P 


pe Ey (A2) 


P 
11 Pay Poo 22 


Equation A2 is the equation of an ellipse in the coordinate 
system (x, y). The ellipse can then be projected on to the 
3 orthogonal planes and the ellipticities, polarization 
angles and sense of polarization on each plane can be 
determined. 

Thus far we have discussed the cross-spectral 
matrices Jaj at a given station for the determination of 
| the polarization parameters. Station-to-station cross- 


spectral estimates Jaj (ri r_, €) (where Jaj is the cross- 


n 
spectral-estimate at frequency f of component i at station 


Bh and component j at station ry) can be calculated. These 
estimates can then be used to determine the coherency and 
phases between pairs of stations according to the following 
expressions. The coherency estimator is given by [see 


Jenkins and Watts, 1968] 


Le: . ti _. is the station-to-station cross- 
where Ji; Li i QG 


spectral estimate. The coherency Kaj is a measure of the 


correlation of component i at station r, with component J at 


station r.. The phase estimator is given by 
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The phase estimator Fay is a measure of the phase difference 
between component i at station re and component j at station 
rae 
The raw auto-power and cross-power spectral esti- 
mates are calculated by means of the Fast Fourier Transform 
algorithm [Gentleman and Sande, 1966] and are subsequently 
smoothed for the determination of the polarization parameters 
and phase estimators. The necessary conditioning of the raw 
data, such as removal of the mean, detrending by passing 
through a digital bandpass filter, cosine-tapering the initial 
and final 10% of the data and addition of zeroes, Ts fist 
carried out before applying the Fast Fourier Transform. 


Let X. be the time series consisting of the original 


N 
time series x, and added zeroes to give N points. The 


discrete Fourier transform of the time series Xn is 


oe i2njn/N n=0,1,25.-.5 


” Xy(5) 
: Ped ; 
where n is a sample at point nA Gee Ge NAt and At is the 
ve 


sample interval). The periodogram is given by 
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I(n) = a | A(n) | 
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The periodogram can be smoothed by taking a simple moving 


average 


‘¢ k . 
VSS i ' E  hntte = par) 


Such smoothing yields a resultant spectral window which has 
roughly a rectangular shape. The shape of such a spectral 
window corresponding to simple moving average smoothing of 
the periodogram for k=7 and m/N = 0.9 is shown in Figure BI: 
The variance of the spectral estimate P depends on the 
bandwidth of the spectral window. An example of the varia- 
tion of the variance with bandwidth is shown in Figure Ba 
It can be seen that one trades variance with resolution and 
vice versa. The spectral estimates used for the computation 
of polarization parameters and phase estimators were smoothed 
using a wide spectral window with 14 degrees of freedom. 


The raw cross-power spectral estimates of two time 


series X(t) and Y(t) are given by 


where A(n) and B(n) are the discrete Fourier transform of 
X(t) and Y(t) respectively. To obtain the smoothed cross- 
power spectral estimates, the real and imaginary parts of 
the spectrum are averaged separately. The spectral window 


in each case is identical to that for the auto-power Case. 
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APPENDIX C 
INFERENCE OF ELECTRIC CURRENTS FROM GROUND MAGNETIC DATA 


Ground-based magnetometers measure the integrated 
effect of all magnetospheric and ionospheric current systems. 
A useful technique to infer the flow of currents from 
ground magnetic data is to portray the magnetic perturba- 
tions in thexform.of,a. latitude profile: ..In the) latitude 
profile, the magnitudes of three components CH Dea Zieare 
plotted as a function of the latitude of the observing 
stations at a given instant of time assuming all the stations 
lie on a common magnetic meridian. Then, by comparing the 
observed latitude profile to theoretically calculated 
latitude profiles based on certain model current systems, 
one can comment on the strength of the current system, the 
borders of the current, etc. This technique has been used 
by Walker [1964], Bonnevier et al. [1970] and Kisabeth and 
Rostoker [1971]. Since our line of stations was positioned 
approximately along a common geomagnetic meridian, displaying 
the magnetic data in the form of latitude profile rendered 
the comparisons between. observed profiles and theoretical 
profiles particularly meaningful. 

Kisabethel 19725 197 bdehasyused 3-dimensional 


equivalent current systems to model the polar electrojet. 
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He formulated Biot-Savart's law in spherical coordinates for 
an arbitrary current loop external to the earth's surface. 
The result of his computations is an expression for a compo- 
nent of the magnetic perturbation (at a point in the surface 
of the earth due to a current loop) in the form of an integral 
of a simple third order matrix which may be evaluated 
numerically for an arbitrary three-dimensional current system 
in the magnetosphere. 

A few samples of theoretical latitude profiles 
computed by Kisabeth [1972] are presented here to illustrate 
how different profiles look for different current systems. 

Consider an ionospheric current system 20°onga and 
5° wide flowing westward with field-aligned Birkeland currents 
flowing down the field lines to the eastern edge and flowing 
up the field lines at the western edge of the. ionospheric 
westward electrojet. Figure Cl shows the latitude profile 
which would be observed on a meridian 6° west of the center 
of the current system. The shaded rectangle in Figure Cl 
represents an ionospheric current of one million amperes 
distributed uniformly across the electrojet. It can be seen 
that the borders of the electrojet are delineated by the 
The central portion of the electrojet is 


extrema in AZ. 


characterized by AZ = 0 and the AH maximum. The D profile 


represents the difference between two large contributions 


from Birkeland currents at the edges of the electrojet. The 
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positive H bay regime and the D-component perturbation 
pattern would be absent had the perturbation been generated 
by a purely westward ionospheric electrojet. 

Figure C2 shows latitude profiles associated with 
current systems having various latitudinal current distri- 
butions. Figure C2A shows a latitude profile which would be 
observed on a meridian 5° east of the center of the iono- 
spheric electrojet having the same width and length as the 
current for Figure Cl. The current density has a normal 
distribution across the electrojet. It can be seen that 
each border of the electrojet lies approximately midway 
between the extremum in AZ and AH = 0. Since this profile 
represents observations east of the center of the electrojet, 
the D profile has signs opposite to that of the D profile in 
Figure Cl which represents observation west of the center of 
the electrojet. Figure C2B shows a profile. 10° east of the 
center of an ionospheric westward electrojet of length 40° 
and width 10° with Birkeland currents at the eastern and 
western edges. The current distribution is enhanced at the 
northern border. The asymmetries in the latitudinal current 
distribution can be observed from the latatude profiles sihe 
effect of the enhanced northern border is exhibited clearly 
in the H and Z profiles. Figure C2C shows the magnetic effect 
xisting eastward and westward electrojets, 


associated with coe 


each having different central meridians. Note the enhanced 
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-Z regime between the two electrojets due to the same sense 

of the Z components of the two electrojets~in that region, 

The eastward electrojet produces a positive H profile opposite 
to that produced by a westward electrojet. Figure C2D shows 

a profile due to a double westward electrojet, each having 

the same parameters as those for the system in Figure C2A. 
Note the partial cancellation of the Z component in the 

region between the two electrojets. 

If a profile is? obtained for a simple three- 
dimensional E-W current system far to the east or west of 
the line of stations, the component profiles are considerably 
broadened and the separation of the peaks in AZ will be 
greater than the width of the actuay current. system. 

Having shown the theoretical latitude profiles due 
to E-Wetonospheric currents with Birkeland currents at the 
eastern and western edges, we now present two theoretical 
profiles computed by Kisabeth [1972] for a N-S three- 
dimensional current system. In this model system the 
Birkeland currents flow from the magnetosphere to the 
northern edge of the southward directed ionospheric sheet 
current and out at the southern edge back to the magnetosphere. 
The N-S current is 4° long and 4° wide. Figure (C3 shows a 
latitude profiterc: east, of the central meridian and Figure C4 
shows a latitude profile 2° west of the central meridian. It 
can be seen that the Z and H profiles are reversed for both 


cases while the D profiles remain the same. 
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There is distinct similarity between the mode | 
latitude profiles of the three-dimensional E-W current 
Systems and observed latitude profiles in the morning sector 
in -the H and Z component. Although the exact locations of 
the borders of the electrojet cannot be determined because 
of lack of knowledge of the exact latitudinal current 
distribution, we can, nevertheless, use the Z component 
extrema to delimit the latitudinal range over which signifi- 
cant Beer o cet current flows. If we define the electrojet 
borders by the extrema in AZ (for uniform latitudinal eur eene 
Heed distribution) or one degree poleward/equatorward of 
extrema in AZ (for normal current density distribution across 
the electrojet), we can estimate the equatorward border of 
the electrojet to within + 1° and the poleward border to 
Wien teles° (the larger probable error for the poleward 
border estimate stemming from the sparser station coverage 
at higher latitudes). 

The model three-dimensional E-W current systems 
used in computing the theoretical profiles described in 
this Appendix restrict the Birkeland current to flow at each 
end of the ionospheric electrojet. In fact Birkeland current 
flow is probably distributed along the entire length of the 
westward (or eastward) electrojet (see Figure 5). Therefore, 
the observed D component profile does not correspond to the 


computed D profile. The observed D profile also provides 
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information regarding the degree of tilt of the current 
system waren respect to the line of observing stations. In 
inet eee the observed D-component profile shows a positive 
step-like function across the electrojet in the morning 
sector [Hughes and Rostoker, 1976]. This is an indication 

of the presence of net downward field-aligned currents noted 
by Yasuhara et al. [1975] using polar orbiter data from the 
TRIAD satellite. 
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